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Today,  the  earth's  instrinsic  magnetic  field  prevents 

the  streaming  solar  wind  plasma  from  directly  interacting 

* 

with  the  terrestrial  atmosphere.  Periodically  the  geomag¬ 
netic  field  reverses  direction,  and  during  the  transition 
period,  the  earth's  magnetic  barrier  is  thought  to  dis¬ 
appear.  In  the  past,  studies  about  the  atmospheric-solar 
wind  interaction  dynamics  were  hindered  by  a  scarcity  of 
observational  data.  However,  since  December  1978,  the 

Pioneer-Venus  Orbiter  has  been  providing  daily  observations 
of  the  atmospheric  dynamics  produced  by  the  direct  solar 
wind  interaction  with  the  atmosphere  of  Venus  (the  only 
planet  known  to  lack  an  intrinsic  magnetic  field).  This 
thesis  develops  the  first  three  dimensional  magnetohydro¬ 
dynamic  ( MHD)  theory  of  this  interaction. 

Within  the  ionosphere  of  Venus  analytic  solutions  to 
the  MHD  equations  are  possible  because  of  a  favorable  geo¬ 
metry  between  the  induced  ionospheric  magnetic  fields  and 


the  ionospheric  plasma  motions.  It  is  shown  that  variations 
in  the  solar  wind  speed  and  interplanetary  magnetic  vector 
direction  cause  variations  in  the  dayside  ionospheric  plasma 
flows  and  the  induced  magnetic  field  configuration,  and  that 
these  changes  can  account  for  the  variety  of  magnetic  struc¬ 
tures  observed  by  Pioneer-Venus.  Portions  of  the  Venus 
ionosphere  are  shown  to  be  susceptible  to  the  Kelvin- 
Helmholtz  shear  instability.  The  unusual  shape  of  the  com¬ 
puted  region  of  stability  is  shown  to  be  an  important  key  to 
understanding  the  highly  variable  Pioneer-Venus  observa¬ 
tions.  Model  calculations  are  compared  to  observations  for 
a  number  of  selected  orbits,  and  the  model  is  shown  to  match 
the  observations  in  fine  detail. 
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ABSTRACT 


A  THREE  DIMENSIONAL  MODEL  OF  THE  PLASMA  FLOW  AND 
MAGNETIC  FIELDS  IN  THE  DAYSIDE  IONOSPHERE  OF  VENUS 


by 


Thomas  F.  Tascione 

Today,  the  earth's  instrinsic  magnetic  field  prevents 
the  streaming  solar  wind  plasma  from  directly  interacting 
with  the  terrestrial  atmosphere.  Periodically  the  geomag¬ 
netic  field  reverses  direction,  and  during  the  transition 
period,  the  earth's  magnetic  barrier  is  thought  to  dis¬ 
appear.  In  the  past,  studies  about  the  atmospheric-solar 
wind  interaction  dynamics  were  hindered  by  a  scarcity  of 
observational  data.  However,  since  December  1978,  the 

Pioneer-Venus  Orbiter  has  been  providing  daily  observations 
of  the  atmospheric  dynamics  produced  by  the  direct  solar 
wind  interaction  with  the  atmosphere  of  Venus  (the  only 
planet  known  to  lack  an  intrinsic  magnetic  field).  This 
thesis  develops  the  first  three  dimensional  magnetohydro¬ 
dynamic  ( MHD )  theory  of  this  interaction. 

Within  the  ionosphere  of  Venus  analytic  solutions  to 
the  MHD  equations  are  possible  because  of  a  favorable  geo¬ 
metry  between  the  induced  ionospheric  magnetic  fields  and 


the  ionospheric  plasma  motions.  It  is  shown  that  variations 
in  the  solar  wind  speed  and  interplanetary  magnetic  vector 
direction  cause  variations  in  the  dayside  ionospheric  plasma 
flows  and  the  induced  magnetic  field  configuration,  and  that 
these  changes  can  account  for  the  variety  of  magnetic  struc¬ 
tures  observed  by  Pioneer-Venus .  Portions  of  the  Venus 
ionosphere  are  shown  to  be  susceptible  to  the  Kelvin- 
Helmholtz  shear  instability.  The  unusual  shape  of  the  com¬ 
puted  region  of  stability  is  shown  to  be  an  important  key  to 
understanding  the  highly  variable  Pioneer-Venus  observa¬ 
tions.  Model  calculations  are  compared  to  observations  for 
a  number  of  selected  orbits,  and  the  model  is  shown  to  match 
the  observations  in  fine  detail. 
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1. 


INTRODUCTION. 


Today,  the  earth's  intrinsic  magnetic  field  prevents 
the  streaming  solar  wind  plasma  from  directly  interacting 
with  terrestrial  atmosphere.  Fossilized  rocks  show  that  the 
geomagnetic  field  reverses  direction  at  intervals  as  short 
as  50,000  years  (Haymes,  1971).  During  the  field  transition 
period,  the  geomagnetic  field  is  thought  to  weaken  enough  to 
allow  the  solar  wind  to  directly  interact  with  the  atmo¬ 
sphere.  In  the  past,  studies  about  the  nature  of  this 
interaction  were  hindered  by  the  lack  of  observational 
data.  However,  evidence  of  such  an  interaction  is  presently 
available  on  Venus,  the  earth's  sister  planet.  Venus  is  the 
only  planet  whose  ionosphere  is  known  to  directly  interact 
with  the  solar  wind.  Since  December  1978,  the  PIONEER  VENUS 
ORBITER  (PV)  has  been  providing  the  most  comprehensive  set 
of  observations  ever  obtained  of  an  extraterrestrial  plan¬ 
etary  ionosphere-atmosphere  system.  Once  every  24  hours, 
the  orbiter  collects  data  within  the  solar  wind  interaction 
region;  the  highly  inclined  orbit  (105°)  requires  about  243 
days  to  precess  through  all  local  times  (Colin,  1979).  The 
purpose  of  this  thesis  is  to  develop  a  three  dimensional 
model  of  the  ionospheric  dynamics  resulting  from  this  solar 
wind/Venus  ionosphere  interaction. 

Since  Venus  lacks  an  intrinsic  magnetic  field  (Russell 
et  al . ,  1980b) ,  the  ionosphere  deflects  the  incident  solar 
wind  by  thermal  pressure  and  by  magnetic  pressure  produced 


2 


by  ionospheric  currents.  In  an  earlier  series  of  papers 
(Cloutier  et  al . ,  1969  ?  Michel,  1971?  Cloutier  and  Daniell, 
1973?  Cloutier  et  al .  ,  1974?  Cloutier,  1976)  it  was  shown 
that  a  planetary  atmosphere  that  is  optically  thick  at 
ionizing  wavelengths  should  deflect  most  of  the  solar  wind 
plasma  around  the  ionosphere  forming  a  detached  bow  shock. 
The  planetary  atmosphere  acts  as  a  hard  obstacle  absorbing, 
at  most,  a  few  percent  of  the  incident  solar  wind  flux 
(Daniell  and  Cloutier,  1977?  Cloutier  and  Daniell,  1979). 
PV  observations  have  confirmed  these  theoretical  predic¬ 
tions. 

Observations  reported  by  Wolfe  et  al .  (1979)  and 

Mihalov  et  al .  (19&0)  confirm  the  existence  of  a  strong 

detached  bow  shock.  The  average  nose  radius  for  the  bow 
shock  is  about  1.38  -  1.46  Rv  (radius  of  Venus)  with  fre¬ 
quent  daily  variations  in  excess  of  .16  Rv  (Slavin  et  al., 
1980?  Theis  et  al.,  1981).  However,  Smirnov  et  al .  (1980) 

report  a  nose  radius  of  1.27  Ry.  There  are  two  reasons  for 
the  reported  range  in  bow  shock  locations.  First,  most 
shock  crossings  are  observed  at  large  solar  zenith  angles 
(SZA),  which  makes  extrapolation  to  the  nose  value  diffi¬ 
cult.  Second,  different  procedures  were  used  to  find 
average  bow  shock  shapes  which  best  fit  the  PV  observations 
(Smirnov  et  al . ,  1980).  The  large  daily  variation  of  the 
bow  shock  is  consistent  with  the  theory  that  only  a  few 
percent  of  the  solar  wind  flux  is  absorbed  (Gombosi  et  al . , 
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1980  and  Mihalov  et  al .  ,  1980).  On  the  downstream  side  of 
the  bowshock  is  the  ionosheath  region  where  the  shocked 
solar  wind  plasma  is  diverted  around  the  ionosphere  of 
Venus.  In  the  ionosheath  region,  the  magnetic  field  drapes 
around  the  planet  and  the  peak  magnetic  field  strength 
occurs  in  a  region  where  the  flow  stagnates  (Elphic  et  al .  , 
1980;  Vaisberg  et  al .  ,  1980).  In  the  stagnation  region  the 
plasma  density  is  low  because  of  plasma  losses  along  the 
magnetic  field  lines  ( Zwan  and  Wolf,  1976). 

The  boundary  between  the  hot,  tenuous  ionosheath 
plasma  and  the  cool,  dense  ionospheric  plasma  is  called  the 
ionopause.  Generally,  the  dayside  ionopause  is  easily  iden¬ 
tified  by  a  relatively  abrupt  cutoff  of  the  thermal  ions 
(Taylor  et  al.,  1979a).  However,  the  ionopause  altitude 
exhibits  large  daily  variations;  the  ionopause  height  is 
usually  confined  below  400  km  near  the  subsolar  point  and 
can  extend  beyond  1000  km  near  the  terminator  (Taylor  e_t 
al . ,  1979b,  1980,  1981).  The  ionopause  boundary  occurs  at 
the  altitude  where  the  ionospheric  thermal  pressure  approxi¬ 
mately  balances  the  magnetic  pressure  just  outside  the  iono¬ 
pause.  Since  the  magnetic  pressure  just  above  the  ionopause 
is  approximately  equal  to  the  solar  wind  dynamic  pressure, 
the  ionopause  height  variability  is  a  consequence  of  the 
variability  of  the  solar  wind  ram  pressure  (Brace  et  al., 
1979a,  1980;  Knudsen  et  al .  ,  1979a;  Taylor  et  al .  ,  1979; 


Elphic  et  al . ,  1980;  Spenner  et  al.,  1980). 


In  the  two  dimensional  ionospheric  model  of  Cloutier 
and  Daniell  (Daniell  and  Cloutier,  1977?  Cloutier  and 
Daniell  1979),  the  solar  wind  interaction  drives  rapid 
horizontal  convection  within  the  ionosphere  of  Venus.  The 
convection  is  driven  by  pressure  gradients  and  electric 
fields  induced  in  the  ionosphere  by  the  flowing  magnetized 
solar  wind  plasma;  these  electric  fields  produce  electric 
currents  within  the  ionosphere  which  connect  across  the 
ionopause  into  the  ionosheath  region.  The  average  volume 
distribution  of  ionospheric  currents  was  determined  by  mini¬ 
mizing  the  volume  Joule  heating  by  a  variational  technique. 
This  model  was  recently  modified  to  include  the  effects:  a) 
plasma  depletion  at  the  stagnation  point,  b)  development  of 
Kelvin-Helmholtz  shear  type  instabilities,  and  c)  variations 
in  the  solar  wind  dynamic  pressure  and  the  interplanetary 
magnetic  field  (Cloutier  et  al . ,  1981).  This  revised  model 
matched  PV  magnetic  field  strength  observations  very  well 
near  the  prime  meridian  plane  and  showed  that  the  observed 
bundles  of  twisted  magnetic  fields  or  "flux  ropes"  (Russell 
et  al . ,  1979a,  1979b;  Elphic  et  al . ,  1980)  were  produced  by 
the  Kelvin-Helmholtz  velocity  shear  instability.  The  model 
also  showed  that  changes  in  the  interplanetary  magnetic 
field  (IMF)  orientation  and  the  solar  wind  dynamic  pressure 
could  have  dramatic  effects  on  the  ionospheric  current  dis¬ 
tribution  which  explained  the  daily  variability  in  the  PV 
magnetometer  measurements.  However,  this  model  was  still 


two  dimensional  and  had  limited  applicability  at  longitudes 
far  from  the  noon  meridian.  Furthermore,  the  model  could 
not  match  the  full  PV  magnetic  vector  observations  in 
detail.  Therefore,  the  next  logical  step  was  to  develop  a 
three  dimensional  model  of  the  dayside  ionospheric  dynamics 
of  Venus. 

The  most  straightforward  three  dimensional  model  is  one 
which  solves  the  complete  set  of  vector  magnetohydrodynamic 
(MHD)  equations.  In  practice  it  is  very  difficult  to  solve 
these  formidable  coupled  second  order  differential  equations 
in  all  but  the  simplest  geometry.  In  fact,  such  a  method 
is  impossible  for  Venus  because  very  few  of  the  necessary 
boundary  conditions  are  known.  Similarly,  this  lack  of 
known  boundary  conditions  makes  a  three-dimensional  varia¬ 
tional  method  impractical  since  there  is  no  way  to  determine 
the  electric  field  variation  with  longitude  without  first 
solving  the  MHD  equations.  Fortunately,  two  pieces  of  ob¬ 
servational  information  help  simplify  the  problem  enough  so 
that  analytic  solutions  for  the  dayside  ionospheric  magnetic 
fields  can  be  found.  The  key  observations  are  that  both  the 
magnetic  fields  (Luhmann  et  al . ,  1980,  1981a)  and  plasma 
flow  (Knudsen  et  al .  ,  1980,  1981)  are  nearly  horizontal  at 
large  solar  zenith  angles.  It  will  be  shown  that  if  both 
the  flow  and  fields  are  horizontal  (but  not  necessarily 
parallel  to  each  other) ,  the  MHD  equations  simplify  and  de¬ 
couple,  and  can  be  solved  in  a  straightforward  manner.  It 
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will  also  be  shown  that  the  plasma  depletion  mechanism  oper¬ 
ating  about  the  stagnation  point  introduces  an  asymmetry  in 
the  flow  geometry,  and  that  this  asymmetry  is  carried  into 
the  ionosphere.  Since  the  ionospheric  conductivity  is  rela¬ 
tively  high,  the  ionospheric  magnetic  fields  are  coupled 
with  the  plasma  flow  and  therefore,  an  asymmetry  also  exists 
in  the  magnetic  field  geometry.  The  above  theory  is  devel¬ 
oped  for  the  dayside  ionosphere  of  Venus  in  Chapter  2;  in 
Chapter  3,  the  horizontal  flow/field  model  calculations  are 
compared  to  PV  observations.  The  thesis  concludes  with  a 
discussion  of  the  model  and  areas  of  possible  model  improve¬ 


ments 


2. 


MODEL  THEORY 


2.1  Horizontal  MHD  Equations 

The  general  equations  describing  the  steady-state, 
three-dimensional  dynamics  of  the  Venus  ionosphere  are: 

V-px  =  Q  -  L  (2.1.1) 

py7Y  =  +  Pa  +  £  X  a  (2.1.2) 

plus  Maxwell's  equations 

7  *  B  =  pc/£0  (2.1.3a) 

V  •  5  =  0  (2.1.3b) 

V  x  B  =  o  (2.1.3c) 

7  *  B  =  u0tl  (2.1.3d) 

where  Q  and  L  are  the  ion  production  and  loss  rates  respect¬ 
ively;  pc  is  the  free  charge  density;  p  is  the  plasma  mass 
density;  y  is  the  plasma  flow  velocity;  g,  g  and  J,  are  the 
electric  field,  magnetic  field  and  current  density;  P  is  the 

thermal  plasma  pressure;  g  is  the  acceleration  of  gravity; 

and  eQ  and  uQ  are  the  dielectric  constant  and  magnetic 
permeability  of  free  space.  The  angle  between  y  and  ft  can 
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be  expressed  as 


v  •  B 

9  =  cos'1  Vv-rtii0  (2.1.4) 

The  following  analysis  implicitly  assumes  a  tangential 
interplanetary  magnetic  field  (IMF). 

The  inertial  coordinate  system  is  shown  in  Figure  2-1. 
The  Z  axis  is  always  aligned  into  the  direction  of  the  solar 

wind  (anti-solar  wind),  the  Y  axis  is  aligned  parallel  to 

the  IMF  direction  and,  therefore,  the  X  axis  is  always  per¬ 
pendicular  to  the  IMF  direction.  The  coordinate  system 
origin  is  at  the  center  of  Venus.  In  terms  of  spherical 

geometry,  m  is  the  polar  angle  and  ^  is  the  azimuthal 

angle.  As  seen  in  Figure  2-1,  surface  distances  are  mea¬ 
sured  in  terms  of  angular  displacements  of  m  and  also 
shown  are  the  surface  unit  vectors.  The  plasma  flow  stream¬ 
lines  are  lines  parallel  to  m.  That  is,  the  model  assumes 
that  the  plasma  flows  radially  (parallel  to  surface)  away 
from  the  plasma  flow  origin. 

The  assumption  that  the  flow  lines  radiate  from  the 
sub-flow  origin  is  based  upon  the  bulk  ion  flow  observations 
of  Knudsen  et  al .  (  1980).  They  report  that,  with  few 
exceptions,  the  flow  is  parallel  to  the  antisolar  direction 
and  generally  increases  in  magnitude  with  altitude.  These 
observations  covered  a  period  of  nearly  two  years.  It  will 


be  shown  that  due  to  the  high  ionospheric  conductivities, 


9 


the  fields  are  closely  coupled  with  the  plasma  flow  and 
therefore,  the  coordinate  system  (tied  to  field)  shown  in 
Figure  2-1  rotates  with  changes  of  the  IMF  orientation. 

Over  a  two  year  period,  the  IMF  varies  regularly  between 

being  oriented  parallel  to  the  ecliptic  plane  to  being 
perpendicular  to  this  plane.  Since  the  orbit  is  fixed 
relative  to  the  ecliptic  plane,  these  IMF  orientation 
changes  mean  that  the  Pioneer-Venus  Orbiter  has  probably 
taken  measurements  at  every  azimuthal  angle,  y.  Therefore, 
the  assumption  of  radial  flowlines  at  all  y  (strictly 
antisolar)  is  a  very  good  first  approximation. 

In  terms  of  these  curvilinear  surface  coordinates,  the 
horizontal  velocity  and  magnetic  fields  are 

A 

y  (r,  m,  y)  =  v  (r,  m,  y)  m  (2.1.5) 

B  (r,  m,  y)  =  Bm(r,  m,  <J0  m  +  B^,  (r,  m,  y)  y  (2.1.6) 

where  both  y  and  §  are  functions  of  height  (r)  and  position 
(m,  y).  The  horizontal  plasma  acceleration  can  be  written 
(Oates,  1974)  as 


Vv  =  X 


3v  - 


7m 


m  - 


2  „ 
-  r 


(2.1.7) 


Substituting  (2.1.5)  into  (2.1.1),  the  continuity  equation 


becomes 
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p  3v  ,  v  3p  _ 
r  Tm  r  '3m 


Q  -  L 


(2.1.8) 


Using  equation  (2.1.7),  the  components  of  the  equation  of 
motion  become 


PV2  -  3P  „„  a.  1  fi  f  B2 


7r 


*  pg  -  4-  ih  (?)  +  ?) 


(2.1.9) 


v  3  v 
p  f  7m 


*  - 1  3P  -  [A  3B*  -  _ 3Bm  +  _ i  ,2  ,  10. 

r  Tin  u  '■r  3m  rsinm  3 ij»  rtanm-*  (2.1.10) 


0  =  - 


1  3  P 

rsinm  Tijj 


B 


_  jm  rl 

ti  L  tr 


3B 


m 


i  3B  , 
1  ♦ 


B 


♦ 


rsinm  3 


rtanm 


(2.1.11) 


where  Ampere's  law  (2.1.3d)  was  used  to  substitute  the  mag¬ 
netic  field  for  the  current  density.  One  additional  equa¬ 
tion  is  obtained  by  substituting  (2.1.6)  into  Gauss's  Law 
(2.1.3b): 


ii-B  +3L_ 
r  3m  m  rtanm 


rsinm  3^ 


0 


(2.1.12) 


Equation  (2.1.8)  through  (2.1.12)  form  the  mathematical 
framework  for  the  horizontal  flow/field  model. 


_M  AG  NET  1C _ 

FIELD  LINE 


f  V 


SOLAR 

WIND 


Figure  2-1.  Coordinate  system.  The  coordinates  are  aligned 
so  that  the  Z  axis  is  in  the  anti-sunward  direction, 
the  X  axis  is  perpendicular  to  the  IMF  direction  and 
the  Y  axis  is  parallel  to  the  IMF.  Also  shown  are  the 
polar  angle  m,  the  azimuthal  angle  f,  and  the  surface 
unit  vectors  m  and  $. 
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2.2  Plasma  Flow  Asymmetry 

It's  easiest  to  analyze  the  plasma  flow  at  the  prime 
meridian  terminator  and  equatorial  plane  terminator  because 
of  published  electron  density  and  temperature  data  near 
these  locations  (Kliore  et  al.,  1979;  Brace  et  al.,  1979; 
Miller  et  al .  ,  1980;  Theis  et  al .  ,  1980).  In  the  following 
analysis,  the  flow  direction  is  referenced  with  respect  to 
the  external  IMF  direction.  That  is,  parallel  and  perpendi¬ 
cular  labels  refer  to  directions  parallel  and  perpendicular 
to  the  IMF  orientation.  In  Figure  2-1,  the  IMF  is  assumed 
to  be  in  the  ecliptic  plane  and,  therefore,  v,  is  within  the 
equatorial  plane  of  Venus  and  is  along  the  noon  meridian. 
We  have  also  assumed  that  within  the  ionosphere,  the  mag¬ 
netic  field  lines  remain  closely  coupled  (long  diffusion 
time — see  Section  2.4)  with  the  plasma  flow.  Figure  2-1 
qualitatively  shows  the  field  lines  being  carried  along  by 
the  ionospheric  plasma  flow. 

In  the  equatorial  plane,  for  horizontal  flow  and 
fields:  v  =  v^  m  and  B  =  Bm  m .  That  is,  within  the  equa¬ 
torial  plane  the  magnetic  field  is  aligned  parallel  to  the 
flow.  Under  steady  state  conditions,  the  continuity  and 
momentum  equations  become: 


3v»  .  v»  ap  n  T 
p  SOT  +  ~  Sm  =  Q  "  L 


J_  =  dP 
r  Tr 


B. 


+  pg  + 


1  (»  (»)  +  JL.) 


(2.2.1) 


(2.2.2) 


13 


v„  av 


p  r- 


II  II 


l  ap 

r  3m 


0  =  - 


ap  l 


(-4-) 


rsinm  a4>  (iorsinm  d<|> 


(2.2.3) 


(2.2.4) 


The  plasma  mass  density  is 


p  =  +  neme  -  nemi 


(2.2.5) 


where  we  have  assumed  the  ionosphere  is  quasi-neutral 
(ne  -  n^).  As  written,  (2.2.5)  is  for  one  ion  species,  but 
it  can  be  easily  generalized  to  any  number  of  ion  species. 
Substituting  (2.2.3)  into  (2.2.1),  we  find 


l  ap  +  v#  dp 
v  r  dm  r  dm 


(2.2.6) 


Equation  (2.2.6)  is  just  a  quadratic  equation  in  v,(  which 
has  the  solution 


(Q 


L)  +  [  (Q 


L)2 


4  ap  apil/2 

r2  ^m  dm-* 

TTp 

r  am 


(2.2.7) 


where  the  gradients  are  evaluated  near  the  parallel  plane 
terminator.  All  the  terms  on  the  right  hand  side  of  (2.2.7) 
are  derivable  from  observations. 

The  velocity  along  the  noon  meridian  (v^  is  derivable 

A  SV 

in  a  similar  way,  with:  B  =  and  14  =  v^m. 


That  is, 
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along  the  perpendicular  direction  (see  Figure  2-1),  B  is 
always  perpendicular  to  v.  The  resulting  equation  is 


(q'-l)  +  [(q  -  l]2  +  i- 


3  P 

3m 


r3P 


B., 


l/2 


v  ,  = 


2 

r 


3  P 

3m 


(2.2.8) 


where  the  gradients  are  evaluated  near  the  perpendicular 
plane  terminator.  The  'least-energy'  calculation  of 

Cloutier  and  Daniell  (1979)  shows  that  along  the  perpendi¬ 
cular  direction,  the  magnetic  gradients  at  the  terminator 
are  very  small  and,  therefore,  they  contribute  very  little 
to  the  magnitude  of'  Vj_.  In  both  equations  (2.2.7)  and 
(2.2.8),  the  plasma  is  assumed  to  be  a  simple  ideal  gas 
where 


P  =  nek  [T0  +  T*)  (2.2.9) 

and  near  the  terminator  Tg  ■  2  T^  (Bauer  et  al . ,  1979;  Brace 
et  al . ,  1979a;  Knudsen  et  al .  ,  1979;  Miller  et  al . ,  1980). 

The  parallel  density  gradients  are  obtained  from  the 
average  electron  density  data  reported  by  Brace  et  al . 
(1979b)  and  Theis  et  al .  (  1980).  Since  the  IMF  is  usually 
near  the  ecliptic  plane  (Ness,  1965;  Brandt,  1970;  Haymes, 
1971;  Hundhausen,  1978),  and  since  the  PV  orbiter  penetrates 
the  Venus  ionosphere  near  the  planetary  equator  (Colin, 
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1980),  then  the  average  observed  electron  densities  should 
be  representative  of  the  density  gradients  parallel  to  the 
IMF.  Indirect  evidence  (Coffey,  1979;  King,  1979;  Mihalov 
et  al . ,  1980;  Slavin  et  al.,  1980)  supports  the  assumption 
that  the  IMF  was  nearly  parallel  to  the  ecliptic  plane  dur¬ 
ing  the  polar  ionosphere  occultation  experiment  of  Kliore  et 
al .  (1979).  Therefore,  the  polar  occultation  measurements 
provide  the  necessary  data  for  estimating  the  density  grad¬ 
ients  along  the  perpendicular  direction.  At  the  terminator, 
the  production  rate  is  approximately  zero  (Q  «*  0)  and  the 
loss  process  is  assumed  to  be  rate  limited  by  charge  ex¬ 
change  reactions  (Nagy  et  al . ,  1979,  1980).  The  analysis  is 
for  a  single  ion  species,  0+,  which  predominates  at  high 
altitudes  (Taylor  et  al.,  1979a,  1979b,  1980)  and  the  re¬ 
action  rates  are  from  Bauer  et  al .  (1979)  and  Nagy  et  al. 
(1980).  Neutral  densities  are  available  from  Nieman  et  al . 
(1979a,  1979b,  1980)  and  von  Zahn  et  al .  (1979,  1980). 
Substitution  of  the  above  information  into  equations  (2.2.7) 
and  (2.2.8)  results  in  the  estimates  of  Vj_  and  v.j  shown  in 
Figure  2-2.  The  range  of  values  for  v±  and  v |(  are  due  to 
chemistry  and  density  gradient  uncertainties. 

The  results  shown  in  Figure  2-2  appear,  at  first,  to  be 
intuitively  incorrect.  One  might  guess  that  flow  perpendi¬ 
cular  to  the  magnetic  field  should  be  restrained  by  field 
tension  and,  therefore,  vL  should  be  smaller  than  v,.  In 
fact,  the  opposite  occurs:  the  magnetic  field  enhances  the 


* 


ujm  ‘3an±inv 


Figure  2-2.  Velocity  profiles.  The  range  of  values  for  the 
parallel  (v,,)  and  perpendicular  (VjJ  velocities  are  due 
to  the  uncertainties  in  the  ionospheric  chemistry  and 
density  measurements.  The  velocity  calculation  is 
arbitrarily  terminated  at  200  km  because  below  this 
altitude  chemistry  dominates  over  dynamics. 


ITY,  km/sec 


perpendicular  flow  velocity.  This  result  is  consistent  with 
the  conclusions  of  the  plasma  depletion  study  of  Zwan  and 
Wolf  (1976).  They  found  that  as  a  flux  tube  passes  through 
the  magnetosheath  (or  ionosheath  for  Venus),  it  experiences 
a  loss  of  plasma  and  an  accompanying  increase  in  magnetic 
pressure  in  order  to  maintain  pressure  balance  with  the 
ambient  environment.  Zwan  and  Wolf  showed  that  the  plasma 
depletion  mechanism  is  not  axisymmetric  and  the  resulting 
velocity  asymmetry  (v,  >  v(  )  increases  in  proportion  to  the 
amount  of  plasma  depletion  (see  Table  2-1). 

-k  it  * 

Table  2-1  lists  the  plasma  depletion  factor  (P  (t  )/p  (0)) 
and  the  velocity  asymmetry  ratio  (vVv;))  as  a  function  of 
the  dimensionless  parameter  t*  which  is  a  measure  of  how 
long  a  flux  tube  resides  near  the  ionopause.  Longer  resi¬ 
dent  times  allow  for  more  plasma  depletion  and,  therefore, 
larger  asymmetries.  Physically,  the  velocity  asymmetry  is 
due  to  the  fact  that  the  complete  pressure  gradient  (gas  and 
fields)  is  available  to  accelerate  the  plasma  perpendicular 
to  the  fields,  whereas  only  the  smaller  gas  gradient  is 
available  to  accelerate  the  plasma  along  the  field  lines. 
That  is,  in  the  depletion  layer,  the  gas  pressure  is  less 
than  (or  equal  to)  the  magnetic  pressure. 

Another  important  feature  of  Figure  2-2  is  that,  within 
the  ionosphere,  the  computed  velocity  asymmetry  is  most  pro¬ 
nounced  at  high  altitudes.  The  maximum  asymmetry  occurs 
around  the  stagnation  region  because  this  is  the  area  where 
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the  highly  asymmetric  ionosheath  plasma  flow  enters  the 
ionosphere.  At  lower  altitudes,  the  velocity  asymmetry  is 
reduced  due  to  frequent  collisions  with  the  ambient,  thermal 
ionospheric  ions.  The  magnitudes  of  the  plasma  flow  veloci¬ 
ties  shown  in  Figure  2-2  are  consistent  with  the  PV  observa¬ 
tions  of  Knudsen  et  al.  (1980,  1981). 

Thus  far,  the  only  assumption  about  the  ionospheric 
plasma  flow  is  that  the  streamlines  are  radial  lines 
centered  at  the  plasma  flow  origin.  Knudsen  et  al.  (1981) 
indicate  that  the  plasma  pressure  gradient  is  the  principal 
accelerating  mechanism  well  below  the  ionopause.  Unfortuna¬ 
tely,  the  lack  of  detailed  azimuthal  density  observations 
makes  it  impossible  to  accurately  compute  the  plasma  accel¬ 
eration  as  a  function  of  solar  zenith  angle.  Therefore,  we 
will  also  assume  that  the  flow  velocity  remains  constant 
along  any  given  streamline.  Physically,  this  assumption 
isn't  too  bad.  Using  the  electron  data  of  Brace  et  al.,  one 
can  compute  the  approximate  pressure  gradients  along  the 
parallel  direction;  by  the  time  the  dayside  plasma  reaches 
the  terminator,  these  pressure  gradients  have  increased  the 
initial  entry  plasma  velocity  (~  2  km/sec)  by  only  about 


20%. 
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TABLE  2-1. +  PLASMA  DEPLETION  FACTOR 


* 

t 

★  *  ★ 

P  (t  )/p  (0) 

vx/v 

0 

1.0 

1.0 

.2 

.84 

1.28 

.41 

.66 

1.72 

.61 

.  48 

2.38 

.82 

.35 

3.22 

1.02 

.27 

4.35 

1.22 

.17 

6.25 

<+  from  R.  A.  Wolf, 


private  communication, 


1980)  . 
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2.3  Magnetic  Field  Equations 

The  horizontal  flow/field  assumption,  combined  with  the 
assumption  that  the  plasma  streams  radially  away  from  the 
sub-flow  region,  produces  the  separation  of  the  magnetic 
field  terms  from  the  plasma  velocity  in  the  equations  of 
motion  (2.1.9-2.1.11).  That  is,  equation  (2.1.11)  used  with 
Gauss's  Law  (2.1.12)  completely  describes  the  magnetic  field 
configuration  and  these  coupled  inhomogeneous  differential 
equations  depend  only  on  the  azimuthal  pressure  gradient  (-^)  . 
The  solution  to  these  differential  equations  is  complicated, 
but  straightforward,  as  shown  in  Appendix  1.  The  resulting 
magnetic  field  equations  are: 


B 

B 


m 

o 


a^  cosm)  sinc|> 
sin2m 


(2.3.1a) 


B 


m 


o 


aQU  -  sini})) 
sinm 


(2.3.1b) 


a  cosm) 
2  ' 


cos<(> 


sin2m 


(2.3.2) 


where  BQ  is  the  maximum  magnetic  field  strength  (stagnation 

point)  and  the  integration  constants  aQ,  a^,  and  a2  will  be 

discussed  later.  The  variable  B  is  due  to  the  effects  of 

m 
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the  azimuthal  pressure  gradient.  Near  the  perpendicular 
direction  (4>  «  0),  the  field  lines  are  nearly  parallel  to 
each  other  (see  Figure  2-1)  and  since  the  plasma  is  free  to 
move  parallel  to  the  magnetic  field,  (which  will  naturally 
minimize  any  azimuthal  pressure  gradient)  then 


lim 
4*  ♦  0 


dP  _ 


0 


(2.3.3) 


We  also  know  that 


lim  3=0 

♦  ♦  0  m  (2.3.4) 


since  along  the  perpendicular  direction,  the  magnetic  field 
is  entirely  parallel  to  £.  Substituting  (2.3.1a)  and 
(2.3.2)  into  (2.1.11)  shows  that 


d<|> 


3_P 

&4» 


(2.3.5) 


which  implies 


lim 

4>  -*■ 


0 


(2.3.6) 


The  form  of  BmQ  given  in  equation  (2.3.1b)  satisfies  both 
the  constraints  of  equation  (2.3.6)  and  Gauss's  Law 
(2.1.12).  A  more  complete  derivation  of  (2.3.1b)  is  given 
in  Appendix  1.  Therefore,  equation  (2.3.1b)  parameterizes 
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the  effects  of  the  azimuthal  pressure  gradient  on  the  mag¬ 
netic  field  geometry  without  requiring  detailed  measurements 
of  these  gradients.  In  fact,  such  detailed  measurements  do 
not  exist  because  the  reported  ionospheric  densities  are 
averaged  by  SZA  without  regards  for  daily  variations  of  the 
magnetic  field  orientation.  From  day  to  day.  Pioneer  Venus 
samples  the  ionosphere  at  various  azimuthal  angles  due  to 
the  IMF  orientation  changes.  However,  information  about  the 
azimuthal  variability  is  lost  in  the  averaging  process. 

Evaluation  of  the  constants  in  equations  (2.3.1)  and 
(2.3.2)  is  complicated  because  all  the  necessary  boundary 
conditions  are  not  known.  Fortunately,  the  constants 
a;  and  a2  can  be  evaluated  by  a  geometric  argument  and  aQ 
can  be  found  from  physical  constraints  imposed  by  observa¬ 
tions.  The  following  discussion  will  only  highlight  the 
arguments  used  to  evaluate  these  constants;  the  detailed 
analysis  is  in  Appendix  2. 

In  Section  2.2,  it  was  shown  that  a  velocity  asymmetry 
exists  within  the  ionosphere  of  Venus:  Vj_  >  v, .  The  sim¬ 
plest  flow  geometry  is  one  in  which  the  spatial  velocity 
distribution  is  elliptical.  That  is, 

v  (r,  m)  v,  (r,  m) 

v  (r,  m,  <10  =  - j-  (2.3.7) 

[vj2  cos2  +  v^2  sin2  it]1'1 

where  v  is  the  flow  speed  at  a  given  azimuthal  angle,  'll. 
Since  along  the  perpendicular  direction  the  ionospheric 
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magnetic  fields  are  closely  coupled  with  plasma  flow  (long 
diffusion  time — see  section  2.4),  then  for  small  i|>,  the 
field  line  shape  must  also  be  elliptical.  At  larger  values 
of  <j»,  the  field  line  shape  changes  because  at  <|»  =  u/2,  the 
ionospheric  field  must  be  parallel  to  the  IMF  direction  (see 
Figure  2-1 ) . 

Earlier,  the  angle  9  was  defined  to  be  the  angle 
between  v  and  B.  Therefore,  at  any  point  along  a  plasma 
streamline,  the  magnetic  field  can  be  defined  as 

B(m,  <|>)  =  B(m,  <j0  cos  0  m  +  3(m,  4»)  sin  9  4>  (2.3.8) 


where  B(m,  4;)  is  the  total  field  strength  at  the  point 
(m,  <)>).  Substituting  (2.3.8)  into  the  third  equation  of 
motion  (2.1.11)  yields 


a_ 

5m 


(in  B ) 


+ 


_ 1_ 

tanm 


1  39 

sinm  3(|> 


(2.3.9) 


where  azimuthal  pressure  gradients  have  been  ignored  be¬ 
cause  <|>  is  assumed  to  be  small  (see  Appendix  1).  The  geo¬ 
metric  constraint  that  the  field  lines  be  nearly  elliptical 
(see  equation  (A. 2.19))  requires  that 

U  =  -  cos  m  ( 1  -  $-)  (2.3.10a) 

bz 


where 
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and 


v±t 

a  =  r  sin  (-p-) 


(2.3.10b) 


b 


r  sin 


(2.3.10c) 


a 

5 


asymmetry  factor 


(2.3. lOd ) 


Substituting  (2. 3.10a — 2.3.10d)  into  (2.3.9)  and  integrating 
gives 


m 

m 


i 


sin  (2  m) 
sin2  (—  m ) 


dm 


(2.3.11) 


Equation  (2.3.11)  is  valid  in  the  limit  that  if)  goes  to  zero 
and  the  integral  is  solvable  for  a  variety  of  asymmetry 
factors  (see  Table  Al-1).  Furthermore,  using  equation 
(2.3.2),  the  left  hand  side  of  (2.3.11)  becomes 


(aj  -  a2  cos  m)  sin2  m^ 
(ai  -  a2  cos  m-jj  sin2  m 


(2.3.12) 


Combining  (2.3.12)  with  (2.3.11),  the  constants  aJ  and  a2 
can  be  solved  relative  to  some  assumed  initial  value  of  B. 
Therefore,  the  two  unknowns,  aL  and  a2,  have  been  replaced 
by  a  single  unknown,  ,  which  is  a  model  input. 

In  order  to  solve  for  the  constant  aQ,  one  needs  to 
determine  how  the  field  strength  varies  with  solar  zenith 
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angle  along  the  parallel  direction.  Vaisberg  et  al .  (1980a, 
1980b)  report  that  the  best  fit  to  the  SZA  variation  of  the 
observed  ionopause  field  strength  is 


Si| ihl  .  COS2  X 


(2.3.13) 


where  BQ  is  the  nose  value  (stagnation  point)  for  B  and  x  is 
the  angle  between  the  ionopause  and  the  solar  wind  direction 
at  a  given  SZA.  However,  equation  (2.3.13)  cannot  be  used 
directly  because  the  PV  observations  used  in  the  Vaisberg 
analysis  were  taken  over  a  period  of  weeks  as  the  orbit  pre- 
cessed  toward  the  terminator.  During  this  time,  there  is  no 
way  to  establish  the  location  of  the  sub-flow  point  and  the 
magnetic  field  orientation.  Therefore,  equation  (2.3.13)  is 
of  limited  value.  This  equation  does  show,  however,  that 
the  field  strength  near  the  terminator  can  be  very  large  and 
it  also  proves  to  be  useful  guide  in  the  determination  of 
the  constant  aQ. 

The  constant  aQ  is  evaluated  by  an  iterative  technique 
described  in  Appendix  2.  The  iteration  method  consists  of 
estimating  an  initial  value  of  B  along  the  parallel  direc¬ 
tion  and  then  computing  (using  and  a2)  both  the  magnetic 
field  lines  and  field  strengths  everywhere  in  the  iono¬ 
sphere.  The  value  of  aQ  is  adjusted  until  the  magnetic 
longitude  variation  of  the  magnetic  field  strength  appro¬ 
ximates  the  behavior  of  equation  (2. 3.) 3).  By  varying  the 
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inputs,  a  library  of  different  ionospheric  flow/field  con¬ 
figurations  is  generated.  Since  the  constant  aQ  represents 
the  effects  of  the  density  gradients  on  the  ionospheric 

field  geometry,  then  by  selecting  various  combinations  of 
aQ ,  a.  and  a2  we  are  effectively  showing  how  the  ionospheric 
magnetic  fields  change  in  response  to  changing  ionospheric 
conditions. 

In  the  proceeding  discussion,  the  terms  'parallel'  and 
'perpendicular'  described  the  ionospheric  field  geometry 
relative  to  the  IMF  orientation.  Given  the  field  geometry 
described  in  this  section,  some  additional  terminology  is  in 
order.  The  parallel  direction  is  the  "magnetic  equator"  of 
the  ionospheric  magnetic  fields  and  the  perpendicular  direc¬ 
tion  is  the  "magnetic  prime  meridian."  This  terminology 
will  be  used  in  the  remainder  of  this  thesis. 


! 
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2.4  Magnetic  Field  Asymmetry 

In  this  section,  we  will  show  that  the  magnetic 
diffusion  time  is  long  compared  to  the  transport  time  for 
the  flow  to  carry  the  field  lines  across  the  surface  of 
Venus.  Therefore,  the  field  lines  are  closely  coupled  with 
the  flow  and  for  small  <l>,  the  magnetic  field  geometry  will 
be  nearly  elliptical.  The  first  step  is  to  determine  the 
diffusion  time  along  the  magnetic  prime  meridian. 

We  define  a  coordinate  system  which  is  unique  to  this 
section  where  the  Z  axis  is  always  parallel  to  the  magnetic 
fields.  In  order  to  simplify  the  analysis,  we  will  assume 
that  the  field  lines  are  nearly  straight  and  that  the  plasma 
pressure  gradients  are  weak.  •  Under  these  conditions,  the 
generalized  Ohm's  Law  (Spitzer,  1963)  becomes 


where 


i  =  3  ‘  (E  +  v  x  b) 


(2.4.1) 


(2.4.2) 


o1  is  the  Pedersen  conductivity,  a2  is  the  Hall  conducti¬ 
vity,  and  <jq  is  the  parallel  or  zero  conductivity.  Invert¬ 
ing  (2.4.1)  and  solving  for  the  electric  field, 
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E  =  a- 1  •  J  -  v  x  B 


(2.4.3) 


where 


,-i  = 


7  x  b 


7  x  [g  1  *  (  1  -  7x  (v  x  b)  = 


(2.4.4) 


Next,  substituting  (2.4.3)  into  Faraday's  Law  gives 

5B 


-  (2.4.5) 


where  Ampere's  Law  (2.1.3d)  was  used  to  substitute  for  the 
current  density,  J.  Since  £  is  aligned  along  the  2  axis, 
equation  (2.4.5)  becomes 


A _ ^ _  [all  +  aiBi  k  .  7  x 

(i0  +  022  ay2  3x2 


(v  X  b) 


M  ft 

at 


(2.4.6) 


For  a  fluid  at  rest,  (2.4.6)  reduces  to  the  diffusion  equa¬ 


tion 
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_i _  ra£§  +  afjBi  _  _  5B 

^o°3  3y2  3x2 J  3t 


(2.4.7) 


where 


a 

3 


(2.4.8) 


and  a3  is  the  Cowling  conductivity.  That  is,  an  initial 
configuration  of  magnetic  field  will  decay  in  a  diffusion 
time 


•c  =  H  a,  L2  (2.4.9) 
o  J 

The  Pedersen  and  Hall  conductivities  for  Venus  were  computed 
by  Daniell  (1976);  these  conductivities  and  the  associated 
diffusion  times  are  listed  in  Table  2-2.  The  diffusion 
distance  (L)  is  6000  km  or  approximately  the  radius  of 
Venus.  Using  a  plasma  flow  velocity  of  approximately 
2  km/sec,  the  travel  time  across  6000  km  is  less  than  one 
hour,  compared  to  a  high  altitude  diffusion  time  of  more 
than  100  hours.  Therefore,  the  magnetic  fields  are  closely 
coupled  with  the  ionospheric  plasma  flow.  Since  the  plasma 
flow  geometry  for  small  <\>  is  elliptical,  then  the  magnetic 
field  geometry  near  the  magnetic  prime  meridian  must  also  be 
approximately  elliptical.  However,  the  field  asymmetry  does 
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not  completely  reflect  the  flow  asymmetry  because  the  field 
asymmetry  does  not  change  with  altitude. 

The  calculations  of  Cloutier  et  al .  (1981)  show  that  at 
high  altitudes,  the  contours  of  constant  magnetic  field 
strength  are  nearly  vertical  lines  along  the  magnetic  prime 
meridian.  This  behavior  is  shown  in  Figure  2-3;  the  field 
strength  contours  are  labeled  as  a  fraction  of  the  total 
field  strength  (i.e.  0.5,  0.4,  etc.).  The  contours  are 
nearly  vertical  because,  as  seen  from  Table  2-2,  the  dif¬ 
fusion  time  decreases  with  increasing  altitude  so  that  the 
field  can  slip  faster  through  the  plasma  at  high  altitudes. 
The  'least  energy'  calculations  of  Cloutier  and  Daniell 
(1979)  show  that  the  slippage  rate  is  such  that  the  field 
strength  remains  relatively  constant  with  altitude,  even- 
though  the  plasma  flow  velocity  increases  with  altitude. 

Along  the  magnetic  equator,  an  exact  calculation  of  the 
scalar  conductivity  requires  detailed  information  about  both 
the  electron  temperature  and  density.  Fortunately,  such  a 
calculation  is  unnecessary  because  the  horizontal  flow-field 
assumption  allows  one  to  solve  the  generalized  Ohm's  Law 
along  the  magnetic  equator.  The  details  of  this  calculation 
are  shown  in  Appendix  3.  The  final  result  of  the  analysis 
(A. 3. 10)  is 


B(h2) 

B(hi )  *  1  (2.4.10) 
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TABLE  2.2.  CONDUCTIVITY  AND  DIFFUSION  TIME 


alt  (km) 

cji  (mho/m) 

a 2 (mho/m) 

200 

.08 

.4 

250 

.  08 

.02 

300 

.01 

IO-3 

350 

7 

x  10"3 

4 

x  IO"4 

400 

6 

x  10-3 

2 

x  10~4 

450 

5 

x  IO-3 

9 

x  IO"5 

500 

4 

x  IO-3 

8 

x  IO"5 

t(sec) 

9.4  x  io7 
3.8  x  io6 

4.6  *  105 

3.2  x  io 5 

2.7  x  105 

2.3  x  IO5 

1.8  x  10  s 


Figure  2-3.  Magnetic  prime  meridian  plane  current  system. 
The  ionopause  shape  follows  that  of  Spreiter  et  al . 
(1970).  Contours  are  labeled  according  to  the  fraction 
of  the  current  lying  outside  that  contour  at  the  local 
longitude.  The  ionosheath  contours  (dashed  lines)  are 
closed  loops  extending  out  to  the  bow  shock.  The 
straight  dashed  lines  withi  '  the  ionosphere  (labeled 
0.1,  1.0,  10.)  are  contours  of  constant  magnetic 

Reynolds  number.  In  the  2-D  model,  the  current 
distribution  is  assumed  to  be  uniform  over  a  wide 
variety  of  longitudes  around  the  prime  meridian  plane, 
and  therefore,  the  contours  also  represent  magnetic 
field  strength  contours.  The  field  lines  are  perpen¬ 

dicular  to  the  plane  of  the  figure  and  the  contour 
labels  also  indicate  the  fraction  of  total  magnetic 
field  strength  along  each  contour. 
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wher'e  the  numerator  is  the  magnetic  field  strength  at  height 
h^,  and  the  denominator  is  evaluated  at  height  h  1  (h2  >  h^; 
the  field  strengths  are  computed  at  the  same  longitude. 
Therefore,  the  contours  of  constant  magnetic  field  strength 
are  vertical  lines  along  both  the  magnetic  equator  and  the 
magnetic  prime  meridian. 

It  is  not  possible  to  evaluate  the  field  strength 
height  variability  at  any  other  azimuthal  angle  and  there¬ 
fore,  we  will  make  the  simplifying  assumption  that,  at  high 
altitudes,  the  ionospheric  field  strength  contours  are  ver¬ 
tical  everywhere.  It  follows  then  that  the  magnetic  field 
asymmetry  factor  must  remain  constant  with  altitude.  Other¬ 
wise,  the  high  altitude  field  lines  would  be  carried  towards 
the  terminator  faster  than  the  low  altitude  field  lines  and, 
under  these  conditions,  the  field  strength  will  not  remain 
constant  with  altitude.  The  exact  field  asymmetry  factor 
will  depend  on  the  physical  conditions  of  the  ionosphere. 
From  day  to  day,  changes  in  the  solar  wind  dynamic  pressure 
and  IMF  orientation  will  cause  changes  in  both  the  iono¬ 
spheric  flow  and  field  configurations.  These  changes  will 
in  turn  effect  the  ionospheric  conductivity  which,  ulti¬ 
mately,  determines  the  amount  of  field  asymmetry.  Further, 
it's  impossible  to  predict  a_  priori  the  field  asymmetry 
factor  for  a  given  day.  In  practice,  the  model  (Chapter  3) 
is  run  for  a  variety  of  field  asymmetry  factors  and  then  the 
best  match  is  found  for  a  specific  set  of  observations. 
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2.5.  Instability  Criterion. 

The  presence  of  a  large  ionospheric  plasma  velocity 
shear  (see  Figure  2-2)  makes  the  Venus  ionosphere  suscepti¬ 
ble  to  magnetohydrodynamic  instabilities,  such  as  the 
Kelvin-Helmholtz  instability.  This  instability  occurs 
(Hesagawa,  1975)  whenever 

k  •  B 

P(k*Y0)2  >  (2.5.1) 

o 

where  p  is  the  plasma  mass  density  which  is  flowing  at  a 
relative  velocity  YQ  with  respect  to  adjacent  plasma  of  the 
same  density,  and  g  is  a  uniform  magnetic  field  extending 
throughout  both  plasma  regions.  Both  g  and  yQ  are  parallel 
to  the  plane  across  which  the  velocity  shear  is  supported. 
Maximum  stability  occurs  for  all  wave  modes  when  YQ  is 
parallel  to  k  and  Q,  and  the  instability  is  most  likely  to 
occur  when  yQ  is  parallel  to  k  and  perpendicular  to  g. 
Equation  (2.5.1)  can  be  simplified  by  assuming  that  all  un¬ 
stable  wave  modes  always  occur  parallel  to  the  flow  direc¬ 
tion  and  therefore, 


o  V  2  >  B*  cos2q 
°  Po 


(2.5.2) 


where  9  is  the  angle  between  yQ  and  g.  As  seen  from  (2.5.1) 
and  (2.5.2),  the  Kelvin-Helmholtz  instability  is  driven  by 
the  kinetic  enrgy  of  the  flowing  plasma  and  is  stabilized  by 
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the  magnetic  tension.  Generalization  of  (2.5.1),  for  a 
finite  thickness  boundary  with  continuous  velocity  shear  and 
vertical  density  gradients,  results  in  an  instability  cri¬ 
terion  similar  to  (2.5.1)  with  the  additional  constraint 
that  the  finite  thickness  boundary  is  stable  for  short  wave¬ 
lengths  (Ong  and  Roderick,  1972).  In  the  flow/field  model, 
the  simple  expression  (2.5.2)  is  used  to  identify  regions  of 
instability,  subject  to  the  wavelength  constraints  of  Ong 
and  Roderick. 

In  an  earlier  paper  (Cloutier  et  al .  ,  1981),  it  was 
shown  that  the  Kelvin-Helmholtz  instability  is  responsible 
for  regions  of  twisted  bundles  of  magnetic  fields,  called 
'flux-ropes,'  which  are  observed  (Russell  et  al . ,  1979a, 
1979b;  Russell  and  Elphic,  1979;  Elphic  et  al . ,  1980)  in  the 
Venus  ionosphere.  A  qualitative  picture  of  the  sequence  of 
events  leading  to  the  formation  of  flux-ropes  is  shown  in 
Figure  2-4  (panels  a-c).  That  is,  a  perturbation  at  the 
shear  boundary  produces  a  wave  which  pushes  the  boundary 
region  into  an  area  of  increasing  plasma  flow  velocity.  The 
velocity  shear  causes  the  wave  crest  to  advance  faster  than 
the  wave  trough,  resulting  in  the  wrapping  up  of  the  field 
lines  into  discrete  bundles  spaced  at  intervals  of  2ir/k.  In 
the  above  example,  yo  is  perpendicular  to  B  which  is  the 
most  unstable  case.  The  bottom  panel  (d)  of  Figure  2-4 


shows  a  schematic  view  of  the  interior  of  a  flux  rope 
(Russell  et  al . ,  1979b)  which  is  representative  of  the  case 
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KELVIN-HELMHOLTZ  INSTABILITY 

rj  OB  ® 


INTERIOR  STRUCTURE  OF  FLUX  ROPE 


(d) 


Figure  2-4.  Sequence  a-c  schematically  shows  the  formation 
of  flux  ropes  by  the  Kelvin-Helmholtz  instability. 
Sample  field  lines  are  labeled  by  the  number  1  to  4. 
Panel  d  shows  the  interior  structure  of  a  flux-rope 
(from  Russell  etal.,  1979b). 
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where  yQ  and  g  are  not  perpendicular.  Since  the  total 
magnetic  flux  within  a  large  volume  containing  the  flux- 
ropes  must  be  conserved,  the  average  field  between  the  flux- 
ropes  must  be  very  small  when  the  field  within  the  flux- 
ropes  is  relatively  large. 

The  actual  amount  of  magnetic  flux  within  the  flux- 
ropes  is  a  function  of  the  azimuthal  angle,  <j>.  The  flux- 
ropes  contain  significant  magnetic  flux  whenever  yo  is 
perpendicular  to  B  since  a  large  number  of  field  lines  are 
bundled  up  in  each  flux  rope.  However,  when  %0  is  parallel 
to  B,  the  magnetic  flux  within  the  flux-ropes  is  zero 
because,  instead  of  gathering  up  field  lines,  the  parallel 
flow  just  rolls  up  the  field  along  the  direction  of  B.  In 
general,  the  PV  orbit  passes  through  regions  which  vary 
between  these  two  extremes  so  that  the  measured  flux-rope 
field  strength  will  often  change  dramatically  along  the 
oribt. 

Cloutier  et  al .  (1981)  also  showed  that  the  flux-rope 
diameter  is  a  function  of  height.  That  is,  the  finite  con¬ 
ductivity  of  the  Venus  ionosphere  is  an  additional  compli¬ 
cation  to  the  formation  of  flux-ropes.  In  order  to  maintain 
flux-rope  formation,  the  convection  velocity  of  the  plasma 
wrapping  up  the  field  must  exceed  the  diffusion  rate  of  the 
field  through  the  plasma.  The  ratio  of  the  plasma  convec¬ 
tion  velocity,  V,  to  the  magnetic  diffusion  rate  is  called 
the  magnetic  Reynolds  number  (Jackson,  1975).  Along  the 
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magnetic  prime  meridian, 

vi 

\,1W-'o,1,l‘  l2-5-3> 

where  a3  is  Cowling  conductivity  (see  Section  2.4)  and  L  is 
the  scale  length  of  the  magnetic  perturbation  (radius  of  the 
flux-rope).  The  straight  dashed  lines  in  Figure  2-3  are 
contours  of  constant  magnetic  Reynolds  number  for  the  prime 
magnetic  meridian  plane;  these  contours  are  labelled  with 
the  value  of  Rm.  Since  <?3  and  are  known  at  each  alti¬ 
tude,  the  contours  of  R^  define  the  scale  size  of  the  flux- 
ropes  formed  at  a  given  altitude.  Further,  since  least- 
action  favors  formation  of  the  smallest  possible  flux-ropes, 
subject  to  the  short  wavelength  constraints  of  Ong  and 
Roderick  (1972),  the  computations  show  that  within  the  un¬ 
stable  regions  the  flux-rope  size  distribution  should  vary 
from  large  flux-ropes  at  high  altitudes  to  small  flux-ropes 
at  low  altitudes. 

A  similar  analysis  can  be  done  along  the  magnetic 
equator  where  the  magnetic  Reynolds  number  now  becomes 


Km,  -  ^o  °o  V»L 
n 


(2.5.4) 


where  the  scalar  conductivity  is 


°o  = 


nee 


mevei  +  + 


(2.5.5) 
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Hanson  (1965)  gives  the  following  empirical  expressions  for 
the  ion  and  electron  collision  frequencies  used  in  equation 
(2.5.5) 


ven  =  5.4  x  10  10  nnTe1/2  (2.5.6) 

»p  3 

vei  =  ^34  +  4,18  log  neTe”3/2  (2.5.7) 

e 

vin=  2.6  x  10-9  (nn+  ni)  M~l/z  (2.5.8) 

where  nn  is  the  neutral  density,  Te  is  the  electron  tempera¬ 
ture,  ne  is  the  electron  density,  n^  is  the  ion  density,  and 
M  is  the  molecular  weight  of  the  ions  and  neutral  particles 
(assumed  to  have  the  same  mass).  Using  the  neutral  densi¬ 
ties  of  Niemann  et  al .  (1980)  and  von  Zahn  et  al .  (1979, 

1980),  and  the  electron  density  and  temperature  data  of 
Theis  et  al .  (  1980  ),  one  finds  that  in  the  dayside  iono¬ 

sphere,  the  electron-ion  collision  frequency  (ve^)  is  more 
than  two  orders  of  magnitude  larger  than  either  the  elec¬ 
tron-neutral  (ven)  or  ion-neutral  collision  (vin)  fre¬ 
quencies.  Therefore,  substituting  (2.5.7)  into  (2.5.5) 


a 

o 


e2T  3/2 


rp  3 

[34  +  4.18  log  [fp— ]  ] 


(2.5.9) 
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Again,  using  the  average  electron  density  and  temperature 
data  of  Theis  et  al .  (1980),  one  finds  that  aQ  varies  slowly 
with  altitude.  Considering  the  uncertainties  of  the  data, 
we  will  assume  that  the  parallel  conductivity  is  constant 
with  altitude.  Since  the  parallel  velocity  increases  with 
altitude,  the  largest  magnetic  Reynolds  numbers  will  occur 
at  high  altitudes  for  large  flux-rope  sizes.  As  discussed 
earlier,  flux-ropes  along  the  magnetic  equator  are  charac¬ 
terized  by  nearly  zero  flux  and  therefore,  observationally , 
their  size  is  poorly  defined. 

In  conclusion,  the  model  would  predict  that  within  un 
stable  regions,  flux-ropes  vary  in  size  from  large  at  high 
altitudes  to  small  at  low  altitudes.  This  prediction  is 
consistent  with  the  flux-rope  distribution  statistics 
reported  by  Russell  et  al .  (1979a,  1979b). 
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2.6  Summary 

The  following  summary  of  equations  form  the  basis  of 
the  computer  model  described  in  Chapter  3.  The  equation 
numbers  are  from  the  section  in  which  the  equations  were 
derived . 

The  equations  describing  the  horizontal  plasma  flow 

are : 


v  =  v  m 


(2.1.5) 


where 


and 


v  = 


vi  va 


(vj  2cos2(i<  +•  v,  2si n2^)1^2 


(2.3.7) 
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(see  Figure  2-2) 
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P  =  nek  (Te  +  Tl)  (2.2.9) 

Q  is  the  ion  production  rate,  L  is  the  ion  loss  rate,  and  0+ 
is  the  predominant  ion  at  high  altitudes.  As  shown  in  Fig¬ 
ure  2-2,  the  velocity  asymmetry  (v^v^  )  increases  with  alti¬ 
tude  . 

The  horizontal  magnetic  fields  are  described  by  the 
following  equations: 


B  =  B_  m  +  B .  4> 


(2.1.6) 


where 


Bm  aQ  (i|>-sini|>)  (a2-a  jcosralsin^ 

BZ  -  s inm  sin^m 


(2.3.1) 


B.  (aj  -  a2cosm) 

=  - ; — 5 -  COS(t) 

Do  sm^m 


(2.3.2) 


where  BQ  is  the  maximum  field  strength  at  the  stagnation 
point.  The  constant  aQ  parameterizes  the  azimuthal  density 
gradient  and  the  constants  a1  and  a2  determine  the  amount  of 
magnetic  field  asymmetry.  The  magnetic  asymmetry  is  con- 
stant  with  altitude. 

Regions  of  instability  occur  whenever 


where  9  is  the  angle  between  v  and  B,  and  is  the  magnetic 
permeability  of  free  space. 

The  above  equations  were  based  on  the  following 
assumptions : 

1)  steady  state  conditions; 

2)  horizontal  plasma  flow  and  magnetic  fields; 

3)  tangential  IMF; 

4)  plasma  flows  radially  away  from  flow  origin; 

5)  flow  speed  is  constant  for  all  SZA; 

6)  unstable  wave  modes  occur  parallel  to  v. 

Taken  together,  assumptions  1,  2  and  4  simplified  the  MHD 
equations  enough  so  that  analytic  solutions  were  possible. 
Assumption  3  simplifies  the  ionosheath,  plasma  depletion  and 
ionospheric  current  analyses.  The  lack  of  detailed  density 
and  pressure  gradient  information  led  to  assumption  5. 
Finally,  assumption  6  greatly  simplifies  the  instability 
analysis. 
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3.  NUMERICAL  MODEL 
3.1  Model  Components 

The  computer  model  consists  of  three  distinct  regions. 
The  largest  region  extends  from  m  (polar  angle)  greater  than 
40°  to  the  terminator,  and  is  called  the  horizontal  region. 
The  other  two  regions  are  the  "sub-flow"  region  (0  <  m  <  10) 
and  the  transition  region  (10  <  m  <  40).  These  regions  are 
shown  in  Figure  3-1  along  with  the  plamsa  flow  streamlines. 

The  horizontal  region  is  characterized  by  both  horizon¬ 
tal  plasma  flows  and  horizontal  magnetic  fields.  The  equa¬ 
tions  describing  the  physics  of  this  region  are  discussed  in 
Section  2.5.  In  this  region,  the  flow  streamlines  are 
radial  lines  centered  at  the  plasma  flow  origin. 

The  "sub-flow"  region  describes  the  area  where  the 
ionosheath  plasma  enters  the  ionosphere.  The  word  "sub¬ 
flow"  is  put  in  quotations  because  this  region  may  or  may 
not  be  aligned  with  the  nose  of  the  bow  shock  where  the 
solar  wind  "flows"  into  the  planetary  obstacle  (Russell  et_ 
al  ♦ ,  1981).  If  the  tangential  IMF  is  orthogonal  to  the 
solar  wind  flow  direction,  then  the  nose  of  the  bow  shock 
and  ionospheric  "sub-flow"  region  are  aligned  (they  also 
coincide  with  the  subsolar  point).  However,  the  ideal  case 
is  rarely  (if  ever)  attained.  Generally,  the  IMF  and  radial 
solar  wind  flow  are  separated  by  some  angle,  <j>  (spiral  angle 
or  garden-hose  angle).  At  Venus,  the  idealized  spiral  angle 
is  about  35a  (see  Figure  3-2).  At  realistic  spiral  angles. 


Figure  3-1.  Model  regions  and  radial  flowlines.  Region  A 
is  the  sub-flow  region  ( 0  <  m  <  10°),  region  B  is  the 
transition  region  (10  <  m  <  40°)  and  region  C  is  the 
horizontal  region  (40  >  m) .  Also  shown  are  the  radial 
flow  lines  which  are  centered  about  the  sub-flow  point; 
the  flow  lines  are  parallel  to  the  planet's  surface  in 
region  C. 
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the  shocked  solar  wind  plasma  has  a  component  of  velocity 
along  the  ionosheath  magnetic  field  lines  which  diverts  the 
sheath  plasma  flow  from  its  initial  radial  (outward  from  the 
Sun)  direction.  The  amount  of  flow  diversion  or  "aberra¬ 
tion"  depends  on  the  solar  wind  flow  speed  and  the  spiral 
angle,  both  of  which  are  highly  variable  from  day  to  day 
(King,  1981).  Therefore,  the  sub-flow  point  will  not 
usually  be  found  at  the  sane  fixed  planetary  location  from 
one  day  to  the  next.  In  the  three  dimensional  model,  we  do 
not  determine  the  exact  location  of  the  sub-flow  point. 
Instead,  the  model  computes  the  center  of  symmetry  of  the 
ionospheric  magnetic  fields  (and  the  flow  pattern),  and  one 
would  expect  that  this  symmetry  center  corresponds  closely 
to  the  physical  sub-flow  point.  A  detailed  calculation  of 
the  sub-flow  location  would  require  not  only  a  sophisticated 
ionosheath  model,  but  would  also  require  data  which  is  not 
available  from  PV  observations. 

In  Figure  3-3,  the  idealized  plasma  flow  streamlines 
are  shown  for  three  different  areas:  1)  outside  the  bow 
shock  (straight  lines);  2)  within  the  ionosheath  (curved 
solid  lines);  and  3)  within  the  ionosphere  (dashed  lines). 
The  shaded  area  corresponds  to  the  "sub-flow"  region  and,  in 
this  region,  the  flow  is  nearly  vertical.  The  weak  hori¬ 
zontal  flow  of  the  sub-flow  region  is  probably  incapable  of 
supporting  the  Kelvin-Helmholtz  type  instability,  and  thus 
the  model  assumes  that  the  "sub-flow"  region  is  always 


49 


stable.  The  magnetic  field  strength  of  the  "sub-flow" 
region  is  determined  from  the  two-dimensional  electrodynamic 
model  (Cloutier  et  al .  ,  1981)  and  is  maintained  at  90%  of 
the  maximum  field  strength  (stagnation  point)  across  the 
entire  "sub-flow"  region. 

The  transition  region  describes  the  area  where  the 
vertically  in-flowing  plasma  changes  to  strictly  horizontal 
flow.  Within  the  transition  region,  the  magnetic  field 
strength  is  again  based  on  the  two-dimensional  electro¬ 
dynamic  model  and  is  given  by 

B(m)  =  B(  40° )  +  AB  *  (1°- j 2  (3.1.1) 


where  AB  is 


jp  =  0.9  -  B(40°)/Bo  (3.1.2) 
Bo 

The  value  of  B(40°)  is  determined  from  the  horizontal  field 
equations  (2.3.1)  and  (2.3.2).  Equations  (3.1.1)  and 
(3.1.2)  are  applied  at  all  azimuthal  angles. 

The  instability  criterion  (2.5.2)  depends  on  the 
angle  9  between  v  and  B.  Within  the  horizontal  flow  region 
(see  equation  2.3.8), 

B1|)(m,<|>)  =  B(m,\|>)  sin  8  (3.1.3) 
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Bm(m,4>)  =  B (m ,  <|> )  cos  0 


(3.1.4) 


and,  therefore, 


LJ  . 

0  =  tan-1  (^1)  (3.1.5) 
m 

That  is,  the  plasma  flows  along  radial  streamlines  (y  =  vm) 
and  on  a  streamline,  the  local  value  of  g  is  given  by  equa¬ 
tions  (3.1.3)  and  (3.1.4)  As  stated  earlier,  the  "sub-flow" 
region  is  assumed  to  be  stable  and,  therefore,  9  is  not 
computed  for  this  region.  In  the  transition  region,  9  is 
approximately  given  by 

9(m)  =  9(40°)  +  (— j-q-- )  (3.1.6) 

Equation  (3.1.6)  is  an  empirical  result  which  provides  the 
necessary  continuity  with  the  results  computed  in  the  hori¬ 
zontal  flow  region. 

One  last  point  about  the  model  concerns  the  shape  and 
definition  of  the  ionopause  surface.  There  have  been  numer¬ 
ous  attempts  (Wolff  et  al.,  1979;  Vaisberg  et  al.,  1980a, b; 
Brace  et  al.,  1980;  Elphic  et  al.,  1980a, b;  Theis  et  al .  , 
1981)  to  understand  the  daily  variability  of  the  dayside 
ionopause  location  in  terms  of  solar  wind  pressure  dynamics. 
These  studies  have  been  moderately  successful  in  describing 
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the  average  ionopause  position,  but  lack  the  ability  to 
describe  daily  ionopause  variations  in  detail.  Therefore, 
the  model  uses  the  in-situ  measurements  by  the  Orbiter  Ion 
Mass  Spectrometer  (Taylor  et  al .  ,  1979a, b,  1980)  to  deter¬ 
mine  the  ionopause  for  each  orbit  analyzed.  Physically, 
these  measurements  are  the  best  way  to  identify  the  thermal 
ion  boundary  (i.e.,  ionopause).  The  OIMS  measurements  pro¬ 
vide  two  measurements  of  the  ionopause  location  for  each 
orbit.  The  remainder  of  the  ionopause  boundary  is  from  a 
"best  fit"  of  the  aerodynamic  ionopause  shape  of  Spreiter 
and  Stahara  (1980)  to  the  two  OIMS  observations.  For  some 
orbits,  the  magnetic  field  gradients  and  ion  gradients  coin¬ 
cide  at  the  ionopause,  while  in  others  they  do  not.  As  will 
be  shown  later,  strong  magnetic  field  gradients  are  not 
reliable  ionopause  indicators. 


3.2  Model  Application — General  Characteristics 

In  this  section  we  identify  symbols  and  parameters 
which  are  common  to  each  set  of  figures  (3  per  set)  used  to 
analyze  the  PV  observations.  The  first  figure  of  the  set 
(e.g.,  Figure  3-5)  shows  the  projection  of  the  model  mag¬ 
netic  field  lines  and  stability  contours  on  the  surface  of 
the  planet.  The  second  figure  (e.g.,  Figure  3-6)  is  an 
ionospheric  cross-section  within  the  plane  of  the  orbit,  and 
finally,  the  third  figure  (e.g.,  Figure  3-7)  contains  a  com¬ 
parison  of  the  model  results  with  observations. 

In  the  first  figure,  the  magnetic  field  lines  are 
represented  by  solid  curved  lines  and  the  field  line  spacing 
is  not  a  measure  of  magnetic  field  strength.  The  region  of 
stability  is  outlined  by  a  set  of  dashed  lines;  the  lowest 
altitude  contour  is  distinguished  by  shading.  The  height  of 
each  contour  is  listed  on  the  bottom  right  hand  corner  of 
the  figure. 

The  typical  shape  of  a  region  of  stability  is  shown  in 
Figure  3-4.  As  expected,  the  main  region  of  stability  sur¬ 
rounds  the  magnetic  equator.  Centered  about  50  degrees  off 
the  magnetic  equator  are  "finger-like"  extensions  of  the 
main  stable  region.  The  shape  and  extent  of  the  stable 
extension  regions  depends  on  the  relative  strength  of  the 
various  terms  making  up  equation  (2.5.2).  During  periods  of 
enhanced  plasma  depletion,  both  the  magnetic  field  strength 
and  magnetic  asymmetry  factor  increase.  The  change  in  the 


asymmetry  factor  reduces  the  angle  9  between  B  and  v. 
Together  these  two  stabilizing  effects  outweigh  the 
destabilizing  effect  of  the  increased  plasma  flow  speed 
produced  by  the  enhanced  plasma  depletion  mechanism. 
Therefore,  periods  of  enhanced  plasma  depletion  are  char¬ 
acterized  by  enlarged  regions  of  stability. 

Another  important  feature  of  the  stable  region  topology 
is  the  rapid  expansion  of  its  horizontal  dimensions  with 
increasing  altitudes.  Once  again  two  competing  factors  are 
at  work:  the  stabilizing  effect  of  decreasing  ion  density 
with  altitude  overcomes  the  destabilizing  effect  of  increas¬ 
ing  plasma  flow  speed  with  altitude.  The  vertical  ion  den¬ 
sity  gradient  is  based  on  the  observations  of  Brace  et  al . 
(1979b),  Taylor  et  al.  (1980),  and  Theis  et  al.  (1980).  The 
vertical  plasma  velocity  gradient  is  taken  from  both  the  ion 
flow  velocity  observations  of  Knudsen  et  al .  (1980)  and  the 
model  velocity  profile  computed  in  Section  2.2.  The  model 
treats  the  mass  density  as  a  constant  for  each  horizontal 
level  below  the  stirong  ionopause  gradients  because  the  sta¬ 
tistical  results  of  Theis  et  al .  (1980)  show  that,  on 
average,  the  electron  density  is  nearly  uniform  over  a  large 
portion  of  the  dayside  ionosphere  (SZA:  0  to  60°).  Theis 
et  al .  (1980)  also  report  little  variation  of  the  electron 
temperature  over  the  same  SZA  range.  Miller  et  al .  (1980) 
report  a  similar  behavior  for  the  ion  densities  and  tempera¬ 
tures  above  200  km.  Therefore,  it  seems  reasonable  to  assume 
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Figure  3-4.  Idealized  stable  reigon.  The  ionopause  is'  at 
350  km  and  the  levels  are  separated  by  50  km.  The 

stable  region  widens  rapidly  with  height  and  is  aligned 
along  the  magnetic  equator  with  its  center  about  the 
plasma  flow  origin.  An  important  feature  is  the 

"notch-like"  unstable  area  between  the  main  stable 
region  and  the  off-axis  stable  extension  region. 
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that,  at  high  altitudes,  there  is  also  little  change  in  the 
ion  flow  velocity  with  SZA  because  of  the  weak  dayside  pres¬ 
sure  gradients. 

In  the  first  figure  of  the  set  (e.g.  Figure  3-5),  the 
center  of  symmetry  for  the  convection  pattern  is  indicated 
by  the  small  solid  circle  and  the  sub-solar  point  is  denoted 
by  the  large  cross.  The  PV  orbit  is  shown  by  the  diagonal 
line;  the  orbit  direction  is  from  the  top  to  the  bottom  of 
the  figure.  The  varying  thickness  of  the  orbit  line  is  used 
to  distinguish  the  time  spent  in  the  ionosheath  (heavy  line) 
from  the  time  the  orbiter  is  within  the  ionosphere.  Once 
inside  the  ionosphere,  orbital  altitudes  (in  kilometers)  are 
shown  by  the  sidebar,  and  the  ionopause  crossings  are  indi¬ 
cated  by  the  arrows  at  the  edges  of  the  sidebar.  The  side- 
bar  may  conceal  a  portion  of  the  field  lines  and/or  stabi¬ 
lity  contours. 

Within  the  ionosphere,  stable  areas  encountered  by  the 
orbiter  are  enclosed  by  a  "box."  If  the  boxed  region  is 
shaded,  it  means  that  the  stable  area  was  encountered  at  (or 
below)  the  level  of  the  lower  (shaded)  stability  contour. 
For  example,  in  Figure  3-5  the  first  stable  area  is  encoun¬ 
tered  shortly  after  the  inbound  ionopause  crossing  (~  350 
km)  .  A  large  unstable  region  extends  from  about  250  km 


inbound  to  about  250  km  outbound;  another  stable  area  is 
encountered  just  before  the  outbound  ionopause  crossing. 
Figure  3-5  also  shows  horizontal  projections  of  the  observed 
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magnetic  field  vectors  which  are  plotted  as  arrows  along  the 
orbit  line.  The  vector  projections  are  shown  only  for 
stable  ionospheric  regions,  and  for  the  entire  ionosheath 
region.  It  is  important  to  remember  that  the  model  is  only 
applicable  within  the  ionosphere  and  therefore,  at  times, 
the  observed  ionosheath  magnetic  field  vectors  may  not  align 
very  well  with  the  model  generated  ionospheric  magnetic 
field  lines. 

The  second  figure  in  the  set  (e.g.  Figure  3-6)  contains 
contours  of  constant  magnetic  field  strength  within  an 
ionospheric  cross-section,  taken  along  the  orbital  plane. 
The  vertical  axis  gives  the  ionospheric  height  (in  kilo¬ 
meters)  while  the  horizontal  axis  shows  the  distance  from 
the  magnetic  equator  (in  degrees).  The  horizontal  scale  is 
centered  at  the  intersection  of  the  orbital  plane  with  the 
magnetic  equator;  positive  distance  corresponds  to  distance 
"above"  the  magnetic  equator  and  negative  distance  is 
"below"  the  equator.  The  ionospheric  field  strength  con¬ 
tours  (solid  lines)  are  relative  to  the  peak  field  strength 
measured  at  the  stagnation  point.  Ionosheath  contours 
(dashed  lines)  are  modifications  of  the  aerodynamic  model  of 
Spreiter  et  al.  (1970)  and  Spreiter  and  Stahara  (1980).  The 
crossmarks  along  the  orbit  show  one  minute  time  intervals 
relative  to  periapsis;  the  orbit  direction  is  shown  by  the 
solid  arrows.  The  sawtooth  pattern  outlines  the  regions  of 
stability  within  the  cross-section  plane.  For  example,  in 
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Figure  3-6,  the  orbit  intersects  the  region  of  stability 
twice.  The  first  intersection  is  with  the  stable  extension 
region  (see  Figure  3-5)  which  appears  as  an  unattached 
sawtooth  pattern  near  30  degrees;  the  orbit  enters  this 
region  on  its  side  and  exits  through  the  bottom.  The  second 
intersection  occurs  as  the  orbit  crosses  the  magnetic  equa¬ 
tor  just  above  250  km.  This  central  stable  region  is  shaded 
below  250  km  in  order  to  coincide  with  stability  contour 
shading  used  in  Figure  3-5. 

As  discussed  in  Section  3.1,  the  ionopause  altitude  is 
determined  from  in-situ  measurements  of  the  Pioneer-Venus 
Orbiter  Ion  Mass  Spectrometer  (OIMS).  The  ion  data  is 
represented  by  the  asterisks  (*)  in  the  top  panel  of  the 
third  figure  (e.g.  Figure  3-7).  Only  the  atomic  oxygen  ion 
data  is  plotted  because  it  predominates  above  200  km  (Taylor 
et  al .  ,  1980,  1981).  The  loss  of  ion  data  near  periapsis 
indicates  that  the  orbiter  has  passed  beneath  the  region  of 
measurable  atomic  oxygen  ion  concentration.  The  arrows 
below  the  time  axis  mark  the  ionopause  locations.  The  top 
panel  also  plots  the  total  observed  magnetic  field  strength 
versus  time;  the  solid  line  represents  measurements  and  the 
dashed  line  represents  the  model.  The  bottom  panel  shows 
how  the  magnetic  field  direction  changes  with  time.  The 
spike-like  field  strength  features  occur  within  unstable 
ionospheric  regions;  they  are  twisted,  narrow  magnetic  field 
filaments  or  "flux  ropes"  (Russell  etal . ,  1979a, b;  Russell 
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and  Elphic,  1979).  The  field  structure  within  flux-ropes 
(Elphic  et  al . ,  1980a;  Wolff  et  al . ,  1980)  produces  the 
rapid  changes  of  field  directions  shown  in  the  bottom 
panel.  Notice  that  the  model  calculations  (dashed  line)  are 
shown  only  for  the  stable  region;  model  results  are  omitted 
for  the  unstable  regions  because  adding  dashed  "spikes" 
would  only  make  the  figure  more  confusing.  For  example,  in 
Figure  3-7,  the  model  predicts  that  a  large  region  of 
instability  extends  from  about  -2.5  minutes  (end  of  dashed 
line)  to  about  2.1  minutes  (start  of  dashed  line). 

One  final  note  about  terminology:  the  phrases  "convec¬ 
tion  center,"  "field  pattern  center,"  "current  system 
center,"  and  "sub-flow  point"  are  all  synonymous  with  the 
phrase  "center  of  symmetry  for  the  convection  pattern." 

In  the  following  sections,  PV  orbits  170,  175,  176, 
177,  188,  and  189  will  be  analyzed.  These  orbits  were 
selected  because  both  the  ion  data  (H.  A.  Taylor,  private 
communication,  1980)  and  full  magnetic  vector  data  (C.  T. 
Russell,  private  communication,  1980)  were  available.  For 
all  of  these  orbits,  the  magnetic  field  deviation  from 
horizontal  is  small  (R.  E.  Daniell,  private  communication, 
1980).  Orbit  175  will  be  discussed  last  because  it  requires 
more  detailed  analysis. 
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3.3  Analysis  of  Orbit  170  (Figures  3-5,  3-6,  and  3-7) 

The  IMF  direction  is  east-southeast  and  the  maximum 
magnetic  field  strength  is  160  gamma.  The  inbound  ionopause 
crossing  is  at  370  km  (-3.75  minutes)  and  the  outbound  cros¬ 
sing  is  at  350  km  (3.6  minutes).  Periapsis  (186  km)  is 
approximately  35°  east  of  the  convection  pattern  center;  the 
convection  pattern  center  is  15°  northwest  of  the  subsolar 
point.  The  magnetic  field  asymmetry  factor  is  2.0. 

Figure  3-5  shows  the  horizontal  projections  of  the 
magnetic  field  and  stability  contours  on  the  surface  of  the 
planet.  The  stability  contours  (dashed  lines)  ares  smaller 
contour  (shaded  region)  outlines  the  250  km  level,  and  the 
larger  contour  delineates  the  350  km  level.  The  orbiter 
enters  an  unstable  area  immediately  after  the  inbound  iono¬ 
pause  crossing.  A  stable  region  is  encountered  at  about  350 
km  and  it  extends  down  to  approximately  275  km  where  the 
ionosphere  again  becomes  unstable.  The  orbiter  remains 
within  an  unstable  ionosphere  until  it  re-enters  the  stable 
region  at  about  270  km  (outbound). 

Figure  3-6  is  a  cross-sectional  view  of  how  the  orbital 
path  intersects  the  region  of  stability.  The  skew  j.r.  the 
field  strength  contours  is  due  to  a  combination  of  geometri¬ 
cal  effects,  the  most  important  being  the  shape  of  the  hori¬ 
zontal  contours  of  constant  magnetic  field  strength  near  the 
magnetic  equator  (see  Figure  4-1).  These  horizontal  con¬ 
tours  are  nearly  parallel  to  the  magnetic  equator  at  large 
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solar  zenith  angles.  Therefore,  the  angle  at  which  the 
orbit  intersects  the  magnetic  equator  is  very  important;  at 
small  intersection  angles  the  orbital  path  will  nearly 
parallel  these  contours,  while  for  large  intersection  angles 
the  orbit  makes  nearly  a  transverse  cut  across  many  of  the 
contours.  Orbit  170  is  somewhere  between  these  two 
extremes. 

In  Figure  3-7,  the  total  magnetic  field  strength  is 
plotted  versus  time.  As  discussed  in  Section  3.2,  the  ion 
data  is  used  to  locate  the  ionopause  crossings  which  are 
marked  by  the  arrows  beneath  the  time  axis.  The  model  field 
strength  (dashed  line)  agrees  well  with  observations.  There 
is  some  discrepancy  near  -3  minutes,  but  disagreement  in 
this  area  is  not  unexpected  because  the  stable  region  occurs 
near  the  point  where  the  field  lines  have  a  significant 
amount  of  curvature.  The  field  line  curvature  rate  is 
determined  by  the  local  ion  gradients  which  are  model 
inputs.  The  field  lines  shown  in  Figure  3-5  were  not  finely 
tuned  for  Orbit  170.  Instead  for  this  orbit — and  for  all 
other  orbits  to  follow — we  used  the  closest  match  from  a 
precomputed  library  which  spans  a  wide  variety  of  iono¬ 
spheric  conditions.  Therefore,  a  better  match  for  Orbit  170 
could  be  found  by  appropriate  adjustments  to  the  ion  gra¬ 
dients.  Since  our  purpose  is  to  show  the  general  applica¬ 
bility  of  the  model,  it  is  not  worthwhile  to  fine  tune  each 
orbit.  As  we  will  see  later,  Orbit  170  is  one  of  the  poorer 
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Figure  3-5.  Orbit  170 — horizontal  projection  of  the  mag¬ 
netic  field  parameters.  The  magnetic  field  lines  are 
represented  by  the  solid  curved  lines.  The  dashed 
lines  represent  the  width  of  the  stable  region  at  two 
altitudes:  the  inner  contour  (shaded)  is  at  250  km  and 

the  outer  contour  is,,  at  350  km  (legend  is  at  bottom  of 
figure).  The  orbital  path  is  shown  by  the  solid  dia¬ 
gonal  line;  the  heavy  line  indicates  that  the  orbit  is 
outside  the  ionosphere.  The  sidebar  is  used  to  show 
the  orbital  altitude  scale  within  the  ionosphere,  and 
the  sidebar  arrows  mark  the  ionopause  crossings.  The 
arrows  along  the  orbit  are  horizontal  projections  of 
the  observed  magnetic  field  vectors.  Within  the  iono¬ 
sphere,  the  arrows  are  shown  only  for  the  stable 
regions  (outlined  by  the  "boxes").  The  solid  circle 
denotes  the  plasma  flow  origin  and  the  large  cross 
represents  the  subsolar  point.  The  scale,  at  the  top 
of  the  figure,  shows  the  radial  distance,  in  degrees, 
from  the  plasma  flow  origin.  The  asymmetry  factor  is 
2.0  and  the  maximum  magnetic  field  strength  is  160 
^ciinind  • 
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Figure  3-6.  Orbit  170 — orbital  plane  cross-section.  The 
ionospheric  contours  of  constant  magnetic  field 
strength  are  shown  by  the  solid  lines;  the  ionosheath 
contours  are  given  by  the  dashed  lines.  The  magnitude 
of  each  contour  is  relative  to  the  maximum  field 
strength  at  the  stagnation  point  (i.e.,  0.5  is  50%  of 
the  maximum) .  The  ionospheric  contours  are  from  the 
three-dimensional  model  and  the  ionosheath  contours  are 
modifications  of  Spreiter's  (1970,  1980)  aerodynamic 

model.  The  sawtooth  pattern  outlines  the  stable 
regions  within  the  cross-section  plane;  the  portion  of 
the  stable  region  below  250  km  is  shaded.  The  cross¬ 
marks  on  the  orbit  are  one  minute  intervals  relative  to 
periapsis.  The  horizontal  scale  is  centered  at  the 
intersection  of  the  orbital  plane  with  the  magnetic 
equator;  positive  distance  (in  degrees)  is  above  the 
magnetic  equator  while  negative  distance  is  below  the 
equator 
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Figure  3-7.  Orbit  170 — comparison  of  the  model  with  obser¬ 
vations.  The  top  panel  shows  the  model  field  strength 
(dash  line)  versus  observations;  model  results  are 
shown  only  for  the  stable  regions.  A  break  in  the  dash 
line  designates  a  region  of  instability  predicted  by 
the  model.  The  field  strength  scale  is  in  gammas.  The 
"spike-like"  field  strength  features  occur  whenever  the 
orbiter  passes  through  twisted  magnetic  filaments 
("flux-ropes")  which  form  within  the  unstable  regions. 
The  asterisks  (*)  are  the  OIMS  measurements  of  the 
atomic  oxygen  ion  densities;  the  density  scale  is  on 
the  far  right  hand  side  of  the  figure.  The  arrows 

beneath  the  time  scale  mark  the  ionopause  crossings. 
The  time  scale  is  relative  to  periapsis.  The  bottom 

figure  shows  how  the  magnetic  field  direction  changes 
with  time;  areas  of  rapid  direction  changes  are  chara¬ 
cteristic  of  the  turbulence  found  in  the  unstable 
regions.  Model  results  are  shown  by  dash  line. 
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matches.  Figure  3-7  also  shows  the  deviation  from  south 
(direction)  for  the  horizontal  component  of  the  magnetic 
field.  As  expected,  the  match  is  poorest  in  the  ionosheath 
region  since  the  model  field  lines  are  valid  only  for  the 
ionosphere . 

The  outbound  ionopause  crossing  for  Orbit  170  is  not 
clearly  defined.  Hartle  et  al .  (1980)  report  that  a  super- 
thermal  layer  of  0+  ions  was  encountered  just  above  300  km 
during  the  outbound  portion  of  the  orbit.  These  superther- 
mal  ions  are  thought  to  be  due  to  a  non-equilibrium  plasma 
region  just  above  the  relatively  unperturbed  ambient  thermal 
ionosphere  (Taylor  et  al.,  1981).  Since  the  ionopause  is 
defined  to  be  the  boundary  of  the  thermal  ion  envelope,  the 
outbound  ionopause  altitude  could  be  lowered  by  about  40  km. 
However,  lowering  the  outbound  ionopause  would  not  change 
the  analysis  because  the  only  effect  would  be  a  reduction  in 
the  size  of  the  large  stable  region  just  below  the  outbound 
ionopause.  That  is,  the  model  would  still  predict  the 
existence  of  stable  magnetic  fields  near  350  km  because  the 
ionosheath  stagnation  region  is  always  assumed  to  be  stable. 
The  ionopause  crossing  time  would  change  from  3.6  minutes  to 
3.3  minutes  (see  Figure  3-7). 
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3.4  Analysis  for  Orbit  176  (Figures  3-8,  3-9,  and  3-10) 

The  IMF  direction  is  nearly  eastward  and  the  magnetic 
field  maximum  is  approximately  190  gamma.  The  inbound 
ionopause  crossing  is  at  280  km  (-3.25  minutes)  and  the 
outbound  crossing  is  at  300  km  (3.5  minutes).  Periapsis 
(145  km)  is  approximately  5°  west  of  the  convection  center; 
the  convection  is  about  20°  northeast  of  the  subsolar  point. 
The  magnetic  field  asymmetry  factor  is  2.5. 

In  Figure  3-8  the  stability  contours  are:  inner  con¬ 
tour  is  for  the  250  km  level,  and  the  outer  contour  outlines 
the  300  km  level.  Orbit  176  remains  within  the  stable 
region  for  the  entire  ionospheric  transit.  A  cross-sec¬ 
tional  view  of  the  orbital  trajectory  through  the  region  of 
stability  is  shown  in  Figure  3-9. 

The  observed  field  strength,  plotted  in  Figure  3-10, 
has  two  areas  of  data  drop  out:  between  1  and  2  minutes, 
and  then  again  between  3  and  4  minutes.  The  PV  fluxgate 
magnetometer  has  an  effective  range  of  ±128  gamma  (Russell, 
1980a).  The  model  field  strength  (dashed  line)  matches  the 
observations  very  well.  The  ion  data  is  again  used  to  iden¬ 
tify  the  ionopause  crossings,  which  are  marked  by  the  arrows 
beneath  the  time  axis.  Figure  3-10  also  shows  a  good  match 
between  the  model  field  direction  and  observations. 
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Figure  3-8.  Orbit  176 — horizontal  projection  of  the  mag¬ 
netic  field  parameters.  The  magnetic  field  lines  are 

shown  by  the  solid  curved  lines.  The  dashed  lines 
represent  the  width  of  the  stable  region  at  two  alti¬ 
tudes:  the  inner  contour  (shaded)  is  at  250  km  and  the 

outer  contour  is  at  300  km  (legend  is  at  bottom  of 

figure).  Within  the  ionosphere,  the  stable  region  is 

enclosed  by  a  "box";  the  shaded  portion  of  the  box 
indicates  that  this  area  is  at,  or  below,  the  250  km 
level.  The  arrows  along  the  orbit  are  horizontal 
projections  of  the  observed  magnetic  field  vectors. 
The  scale  in  the  upper  right  hand  corner  of  the  figure 
is  the  radial  distance,  in  degrees,  from  the  plasma 
flow  (solid  circle).  The  cross  identifies  the  subsolar 
point.  The  asymmetry  factor  is  2.5  and  the  maximum 
field  strength  is  190  gamma. 
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Figure  3-9.  Orbit  176 — orbital  plane  cross-section.  The 
ionospheric  contours  of  constant  magnetic  field 
strength  (solid  lines)  are  from  the  three-dimensional 
model  and  the  ionosheath  contours  (dash  lines)  are 
modifications  of  Spreiter's  aerodynamic  model.  The 
field  strength  magnitudes  are  relative  'to  the  stagna¬ 
tion  point  maximum.  The  sawtooth  pattern  outlines  the 
stable  region  within  the  orbital  plane;  the  portion  of 
the  stable  region  below  250  km  is  shaded.  The  cross¬ 
marks  on  the  orbit  are  one  minute  intervals  relative  to 
periapsis.  The  horizontal  scale  is  centered  at  the 
intersection  of  the  orbital  plane  with  the  magnetic 
equator;  positive  distance  (in  degrees)  is  above  the 
magnetic  equator  and  negative  distance  is  below  the 
equator. 


Figure  3-10.  Orbit  176 — comparison  of  model  with  obser¬ 
vations.  The  top  figure  compares  model  field  strength 
(dash  line)  with  observations;  the  field  strength  scale 
is  in  gammas.  The  aterisks  (*)  are  the  OIMS  measure¬ 
ments  of  the  atomic  oxygen  ion  densities;  the  density 
scale  is  on  the  far  right  hand  side  of  the  figure.  The 
arrows  beneath  the  time  scale  mark  the  ionopause  cross¬ 
ings.  Time  is  relative  to  periapsis:  negative  time  is 
pre-periapsis  and  positive  time  is  post-periapsis.  The 
bottom  figure  shows  the  magnetic  field  direction 
changes  with  time;  model  results  are  given  by  dash 
line . 
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3.5  Analysis  of  Orbit  177  (Figures  3-11,  3-12,  and  3-13) 

The  IMF  direction  is  northeast  and  the  maximum  magnetic 
field  strength  is  110  gamma.  The  inbound  ionopause  crossing 
is  at  350  km  (-4  minutes)  and  the  outbound  crossing  is  also 
at  350  km  (4  minutes).  Periapsis  (150  km)  is  about  5°  west 
of  the  convection  pattern  center  and  the  convection  center 
is  approximately  20°  northeast  of  the  subsolar  point.  The 
magnetic  field  asymmetry  factor  is  1.5. 

As  seen  in  Figure  3-11,  the  orbiter  enters  the  stable 
region  immediately  after  the  inbound  ionopause  crossing;  the 
stable  region  extends  for  about  100  km  (~  1  minute).  The 
orbiter  then  enters  an  unstable  area  which  extends  down  to 
about  an  altitude  of  200  km,  where  again  the  orbiter  enters 
the  region  of  stability.  This  pattern  is  closely  duplicated 
during  the  outbound  portion  of  the  orbit. 

The  topology  of  the  stable  region  is  more  apparent  in 
Figure  3-12.  As  discussed  in  Section  3.2,  the  stable  region 
is  characterized  by  rapid  broadening  at  high  altitudes  and 
it  is  this  "ledge-like"  protrusion  that  Orbit  177  encounters 
at  both  ionopause  crossings.  Once  again.  Figure  3-13  shows 
that  the  model  agrees  well  with  both  the  measurements  of 
field  magnitude  and  direction  within  the  ionosphere.  The 
ionopause  crossings  are  indicated  by  the  arrows  beneath  the 
time  axis. 


Figure  3-11.  Orbit  177 — horizontal  projection  of  magnetic 
field  parameters.  The  magnetic  field  lines  are  shown 
by  the  solid  curved  lines.  The  dash  lines  represent 
the  width  of  the  stable  region  at  two  altitudes:  the 
inner  contour  (shaded)  is  at  250  km  and  the  outer  con¬ 
tour  is  at  350  km  (legend  at  bottom  of  figure).  The 
sidebar  shows  the  orbital  altitude  within  the  iono¬ 
sphere;  the  sidebar  arrows  mark  the  ionopause  cross¬ 
ings.  Within  the  ionosphere,  stable  regions  are 
enclosed  by  a  "box";  the  shaded  portion  of  the  box 
indicates  that  this  area  is  at  or  below  the  250  km 
level.  The  arrows  along  the  orbit  line  are  the 

observed  magnetic  field  vector.  Within  the  ionosphere, 
the  vectors  only  appear  in  the  stable  regions.  The 
scale  in  the  upper  right  hand  corner  of  the  figure  is 
the  radial  distance,  in  degrees,  from  the  plasma  flow 
origin  (solid  circles).  The  cross  designates  the  sub¬ 
solar  point.  The  asymmetry  factor  is  1.5,  and  the 

maximum  field  strength  is  110  gamma. 
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Figure  3-12.  Orbit  177 — orbital  plane  cross-section.  The 
ionospheric  contour  of  constant  magnetic  field  strength 
(solid  lines)  are  from  the  three-dimensional  model  and 
the  ionosheath  contours  (dash  line)  are  modifications 
of  Spreiter's  aerodynamic  model.  The  field  strength 

magnitudes  are  relative  to  the  stagnation  point  maxi¬ 
mum.  The  sawtooth  pattern  outlines  the  stable  region 
within  the  orbital  plane?  the  shaded  portion  is  below 
250  km.  The  horizontal  scale  is  centered  at  the  inter¬ 
section  of  the  orbital  plane  with  the  magnetic 
equator.  Positive  distance  (in  degrees)  is  above  the 
magnetic  equator,  and  negative  distance  is  below  the 
equator. 
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Figure  3-13.  Orbit  177 — comparison  of  model  with  observa¬ 
tion.  The  top  figure  shows  the  model  field  strength 
(dash  line)  versus  observations;  model  results  are 
shown  only  for  the  stable  regions.  A  break  in  the  dash 
line  designates  a  region  of  instability  predicted  by 
the  model.  The  field  strength  scale  is  in  gammas.  The 
"spike-like"  field  strength  features  occur  whenever  the 
orbiter  passes  through  twisted  magnetic  field  filaments 
("flux-ropes")  which  form  in  the  unstable  regions.  The 
asterisks  (*)  are  the  OIMS  measurements  of  the  atomic 
oxygen  ion  densities.  The  arrows  beneath  the  time 
scale  mark  the  ionopause  crossings.  The  bottom  figure 
shows  how  the  magnetic  field  direction  changes  with 
time;  areas  of  rapid  direction  changes  are  character¬ 
istic  of  the  turbulence  found  in  the  unstable  regions. 
Model  results  are  shown  by  the  dash  line. 
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3.6  Analysis  for  Orbit  188  (Figures  3-14,  3-15,  and  3-16) 

This  orbit  is  particularly  interesting  because  the 
orbiter  encounters  an  entirely  unstable  ionosphere.  The  IMF 
direction  is  southeast  and  the  maximum  magnetic  field 
strength  is  110  gamma.  The  inbound  ionopause  crossing  is  at 
370  km  (-4  minutes)  and  the  outbound  crossing  is  at  280  km 
(3  minutes).  Periapsis  (168  km)  is  approximately  26°  north¬ 
east  of  the  convection  pattern  center  and  the  convection 
center  is  about  14°  southwest  of  the  subsolar  point.  The 
magnetic  field  asymmetry  factor  is  1.5. 

As  seen  in  Figure  3-14,  the  orbiter  remains  in  an  un¬ 
stable  portion  of  the  ionosphere  during  the  entire  iono¬ 
spheric  passage;  the  outbound  ionopause  crossing  occurs 
before  the  orbiter  enters  the  stable  region  which  encircles 
the  convection  center.  The  relative  position  of  the  orbital 
path  versus  the  region  of  stability  is  more  apparent  in 
Figure  3-15.  As  in  all  the  previous  figures,  the  OIMS  mea¬ 
surements  are  used  to  determine  the  ionopause  locations. 
However,  the  inbound  ionopause  crossing  for  Orbit  188  is  not 
clearly  defined  and  could  be  higher  than  the  370  km  used  in 
the  analysis.  The  uncertainty  is  due  to  a  high  altitude 
superthermal  component  of  the  0+  ions. 

As  seen  in  Figure  3-16,  the  model  field  strength  and 
direction  agree  very  well  with  observations.  The  model 
predicts  that  the  unstable  region  extends  from  approximately 
-3.5  minutes  to  about  2.5  minutes,  which  again  agrees  with 


observations.  In  the  lower  panel  of  Figure  3-16  there 
appears  to  be  a  stable  region  near  -2.5  minutes  {i.e., 
steady  field  direction)  which  would  disagree  with  the  model 
predictions.  However,  the  field  direction  information  is 
only  part  of  the  available  data;  the  upper  panel  of  Figure 
3-16  shows  that  at  -2.5  minutes  the  field  strength  signature 
is  characteristic  of  an  unstable  region.  This  supports  the 
model's  contention  that  during  this  time,  the  orbit  was 
within  a  region  of  instability. 
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Figure  3-14.  Orbit  188 — horizontal  projection  of  the  mag¬ 
netic  field  parameters.  The  magnetic  field  lines  are 
shown  by  the  solid  curved  lines.  The  dash  lines  repre¬ 
sent  the  width  of  the  stable  region  at  two  altitudes: 
the  inner  contour  (shaded)  is  at  250  km  and  the  outer 
contour  is  at  300  km  (legend  is  given  at  bottom  of 
figure).  The  sidebar  shows  the  orbital  altitude  within 
the  ionosphere;  the  sidebar  arrows  mark  the  ionopause 
crossings.  The  arrows  along  the  orbit  line  are  the 
observed  magnetic  field  vectors.  Since  the  ionosphere 
is  totally  unstable,  these  vectors  are  shown  only  for 
the  ionosheath  (heavy  line).  The  scale  in  the  upper 

right  hand  corner  of  the  figure  is  the  radial  distance, 
in  degrees,  from  the  plasma  flow  origin  (solid  cir¬ 
cles).  The  cross  designates  the  subsolar  point.  The 

asymmetry  factor  is  1.5,  and  the  stagnation  point  field 
strength  is  110  gamma. 
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Figure  3-15.  Orbit  188 — orbital  plane  cross-section.  The 
ionospheric  contours  of  constant  magnetic  field 
strength  (solid  lines)  are  from  the  three-dimensional 
model  and  the  ionosheath  contours  (dash  line)  are 
modifications  of  Spreiter’s  aerodynamic  model.  The 
field  strength  magnitudes  are  relative  to  the  stagna¬ 
tion  point  maximum.  The  sawtooth  pattern  outlines  the 
stable  region  within  the  orbital  plane;  the  shaded 
portion  is  below  250  km.  The  horizontal  scale  is 
centered  at  the  intersection  of  the  orbital  plane  with 
the  magnetic  equator. 
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Figure  3-16.  Orbit  188 — comparison  of  model  with  observa¬ 
tions.  The  top  figure  shows  the  modeled  field  strength 
(dash  line)  versus  observations;  model  results  are 
shown  only  for  the  stable  regions.  A  break  in  the  dash 
line  designates  a  region  of  instability  predicted  by 
the  model.  The  field  strength  scale  is  in  gammas.  The 
"spike-like"  field  strength  features  occur  whenever  the 
orbiter  passes  through  twisted  magnetic  field  filaments 
("flux-ropes")  which  form  in  the  unstable  regions.  The 
asterisks  (*)  are  the  OIMS  measurements  of  the  atomic 
oxygen  ion  densities.  The  arrows  beneath  the  time 
scale  mark  the  ionopause  crossings.  The  bottom  figure 
shows  how  the  magnetic  field  direction  changes  with 
time;  rapid  direction  changes  are  characteristic  of  the 
turbulence  found  in  the  unstable  regions.  The  model 
results  are  shown  by  the  dash  line. 


3.7  Analysis  for  Orbit  189  (Figures  3-17,  3-18,'  and  3-19) 

The  IMF  direction  is  nearly  northward  and  the  maximum 
field  strength  is  140  gamma.  The  inbound  ionopause  crossing 
is  at  550  km  (-5.7  minutes)  and  the  outbound  ionopause 
crossing  is  at  330  km  (3.5  minutes).  Periapsis  (170  km)  is 
approximately  30°  northeast  of  the  convection  pattern  center 
and  the  convection  center  is  10°  southwest  of  the  subsolar 
point.  The  asymmetry  factor  is  2.0. 

The  relatively  high  inbound  ionopause  crossing  and  the 
lower  outbound  ionopause  crossing  is  consistent  with  the 
orbiter  entering  the  ionosphere  relatively  far  from  the 
convection  center  and  exiting  near  the  convection  pattern 
center.  In  the  model,  the  inbound  ionopause  crossing  is 
about  32°  from  the  convection  center,  and  the  outbound 
ionopause  crossing  is  within  10°  of  the  convection  center. 
The  orbiter  enters  an  unstable  area  almost  immediately  after 
the  inbound  ionopause  crossing  and  remains  within  the  un¬ 
stable  area  until  it  just  grazes  the  edge  of  the  main  region 
of  stability  at  about  -4  minutes  (370  km).  There  is  some 
uncertainty  about  how  long  the  orbiter  remains  within  the 
stable  region  because,  at  high  altitudes,  the  model  is  very 
sensitive  to  small  changes  in  the  orbit  orientation  relative 
to  the  stable  region  geometry.  That  is,  the  horizontal  width 
of  the  stable  region  changes  dramatically  with  altitude  (see 
Section  3.2)  and  therefore,  for  grazing  intersection,  small 
shifts  in  the  orbit's  position  relative  to  the  region  of 
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stability  can  drastically  change  the  transit  time  through 
the  region  of  stability.  In  Figure  3-17,  the  field  (and 
stable  region)  is  oriented  so  as  to  give  the  best  fit  with 
observations.  Within  the  orbital  cross-section  plane  (see 
Figure  3-18),  this  stable  region  would  appear  as  a  thin, 
irregular  ellipsoid  bordering  the  orbit  line  from  approxi¬ 
mately  -4  minutes  to  about  -3.5  minutes.  In  the  scale  of 
Figure  3-18,  this  small  ellipsoid  would  be  slightly  wider 
than  the  orbit  line  and  therefore,  it  is  omitted  for  the 
sake  of  clarity. 

As  seen  in  Figure  3-19,  the  model  overestimates  the 
field  strength  at  the  outbound  ionopause  crossing.  Con¬ 
sidering  that  the  ionosheath  model  {Spreiter  et  al .  ,  1970; 
Spreiter  and  Stahara,  1980)  was  simply  grafted  onto  the 
ionospheric  model,  it  is  not  surprising  that,  occasionally, 
spurious  magnetic  field  strength  magnitudes  occur  near  the 
ionopause  boundary.  The  biggest  field  direction  discrepancy 
occurs  at  -4  minutes  where  the  field  lines  have  a  large  cur¬ 
vature.  As  discussed  earlier  (Section  3.3),  "fine-tuning" 


model  inputs  could  improve  the  match. 
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Figure  3-17.  Orbit  189 — horizontal  projection  of  the  mag¬ 
netic  field  parameters.  The  magnetic  field  lines  are 
shown  by  the  solid  curved  lines.  The  dash  lines  repre¬ 
sent  the  width  of  the  stable  regions  at  two  altitudes: 
the  inner  contour  (shaded)  is  at  250  km  and  the  outer 
contour  is  at  450  km  (legend  is  at  bottom  of  figure). 
The  sidebar  shows  the  orbital  altitude  within  the  iono¬ 
sphere;  the  sidebar  arrows  mark  the  ionopause  cross¬ 
ings.  The  stable  region  within  the  ionosphere  is 

enclosed  by  a  "box";  the  stable  region  starts  at  about 
370  km  and  extends  down  to  about  320  m.  The  arrows 

along  the  orbit  line  are  the  observed  magnetic  field 
vectors;  within  the  ionosphere,  the  vectors  appear  only 
in  the  stable  region.  The  scale  in  the  upper  right 
hand  corner  of  the  figure  is  the  radial  distance,  in 
degrees,  from  the  plasma  flow  origin  (solid  circles). 
The  crosshatch  area,  in  the  upper  right  hand  corner, 
designates  the  boundary  of  the  planetary  disk.  The 
cross  denotes  the  subsolar  point.  The  asymmetry  factor 
is  2.0,  and  the  maximum  field  strength  is  140  gamma. 
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Figure  3-18.  Orbit  189 — orbital  plane  cross-section.  The 
contour  of  constant  magnetic  field  strength  within  the 
ionosphere  are  shown  by  the  solid  lines;  ionosheath 
contours  are  given  by  dash  lines.  The  field  strength 
magnitudes  are  relative  to  the  stagnation  point  maxi¬ 
mum.  The  sawtooth  pattern  outlines  the  stable  region 
surrounding  the  plasma  flow  origin.  The  stable  region 
encountered  during  inbound  is  omitted  for  simplicity 
(see  text).  The  horizontal  scale  is  centered  at  the 
intersection  of  the  orbital  plane  with  the  magnetic 
equator. 
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Figure  3-19.  Orbit  189 — comparison  of  model  with  observa¬ 
tion.  Model  results  (dash  line)  are  shown  only  for  the 
stable  regions.  A  break  in  the  dash  line  designates  a 
region  of  instability  predicted  by  the  model.  The 
field  strength  scale  is  in  gammas.  The  "spike-like" 
field  strength  features  occur  whenever  the  orbiter 
passes  through  twisted  magnetic  field  filaments  ("flux- 
ropes")  which  form  in  the  unstable  regions.  The  aster¬ 
isks  (*)  are  the  OIMS  measurements  of  the  atomic  oxygen 
ion  densities.  The  arrows  beneath  the  time  scale  mark 
the  ionopause  crossings.  Regions  of  rapid  direction 
changes  are  characteristic  of  the  turbulence  found  in 
the  unstable  regions. 
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3.8  Analysis  for  Orbit  175  (Figures  3-20,  3-21,  3-22,  3-23, 

and  3-24) 

During  this  orbit,  the  IMF  appears  to  be  undergoing  a 
slow  rotation.  This  interpretation  is  based  on  the  follow¬ 
ing  facts.  First,  Orbit  175  is  the  only  orbit  where  the  IMF 
direction  after  the  ionospheric  exit  does  not  return  to  the 
pre-ionospheric  transit  conditions  (compare  Figure  3-24  with 
Figures  3-7,  3-10,  3-13,  3-16,  and  3-19).  Secondly,  the 
inbound  ionopause  crossing  occurs  nearly  2  minutes  before 
the  appearance  of  strong  magnetic  field  gradients.  This  in 
itself  is  not  unusual,  but  taken  with  the  observed  field 
direction  and  the  location  of  the  region  of  stability,  it 
becomes  impossible  to  match  the  observations  with  a  single 
static  model  representation.  An  IMF  rotation  rate  of  about 
2.5°  per  minute  is  required  for  the  model  to  match  the 
observations;  this  rotation  rate  is  compatible  with  in-situ 
observations  (Brandt,  1970;  Hundhausen,  1972).  It  may  be 
that  the  rotation  is  unnecessary  because  the  model  might  not 
be  applicable  for  Orbit  175.  However,  considering  the 
model's  success  with  all  the  other  orbits,  we  will  pursue 
the  idea  that  the  IMF  is  rotating. 

The  model  assumes  that  both  the  ionospheric  magnetic 
fields  and  plasma  flow  are  strictly  horizontal.  Observa¬ 
tions  indicate  that  for  Orbit  175,  the  magnetic  field 
remained  horizontal  during  the  entire  orbiter  transit. 
However,  the  rotating  magnetic  field  will  induce  an  electric 
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field  which  can  accelerate  ions  to  speeds  comparable  to  the 
ambient  flow  speed,  a  non-ignorable  effect.  If  "frozen- in" 
flux  conditions  exist  in  the  ionosheath  region,  then  this 
additional  plasma  acceleration  increases  the  plasma  flow 
speed  perpendicular  to  the  local  ionosheath  magnetic  field. 
The  net  affect  is  an  increase  in  both  the  plasma  flow  and 
magnetic  field  asymmetry  factors.  The  analysis  assumes  that 
the  plasma  depletion  mechanism  is  unaffected  by  the  magnetic 
field  rotation. 

The  inbound  ionopause  crossing  is  at  about  450  km  (-4.8 
min)  and  the  outbound  ionopause  crossing  is  at  approximately 
300  km  (3.5  min).  Periapsis  (148  km)  is  nearly  25°  north  of 
the  convection  pattern  center  and  the  convection  center  is 
about  20°  east-southeast  of  the  subsolar  point.  The  convec¬ 
tion  center  remains  fixed  during  the  rotation.  The  maximum 
field  strength  is  160  gamma  and  the  magnetic  field  asymmetry 
factor  is  2.0. 

Figures  3-20,  3-21,  and  3-22  show  the  model  analysis  at 
6  minute  intervals  during  the  rotation.  Each  figure  is  an 
average  representation  for  the  altitude  range  shown  in 
brackets  on  the  right  hand  side  of  the  figure.  For  example. 
Figure  3-20  is  an  average  representation  for  the  altitude 
range  of  450  to  300  km.  As  usual,  a  sidebar  is  used  to  show 
the  altitude  scale  within  the  ionosphere.  Figure  3-21 


covers  the  altitude  range  of  250  to  148  km;  no  magnetic 
field  vectors  are  shown  because  the  ionosphere  is  unstable 


over  this  altitude  range.  The  outbound  magnetic  field  con¬ 
figuration  is  shown  in  Figure  3-22.  In  this  figure,  the 
orbit  remains  within  the  region  of  stability  except  for  a 
short  time  when  it  crosses  the  narrow  unstable  area  between 
the  main  stable  region  and  the  stable  extension  region.  The 
ionosheath  magnetic  field  vectors  are  shown  for  complete¬ 
ness. 

Figure  3-23  is  a  composite  picture  of  the  orbital 
cross-section  plane.  The  symmetry  of  the  magnetic  field 
strength  contours  ("above"  and  "below"  the  magnetic  equator) 
is  an  artifact  of  the  averaging  process.  The  stable  regions 
are  not  shown  in  this  figure  because  it  is  impossible  to 
depict  the  relative  separation  of  these  stable  regions  in  a 
composite  format.  Once  again,  the  ion  data  is  used  to 
determine  the  ionopause  altitude.  As  seen  in  Figure  3-24, 
the  inbound  ionopause  crossing  is  not  near  the  observed 
strong  magnetic  field  gradient.  Instead,  the  strong  mag¬ 
netic  field  gradient  coincides  with  the  boundary  between  the 
stable  and  unstable  regions.  Of  course,  if  the  ionosphere 
is  unstable  immediately  inside  the  ionopause,  then  the  iono¬ 
pause  crossing  will  be  characterized  by  a  strong  magnetic 
field  gradient.  The  important  point  is  that  strong  magnetic 
field  gradients  are  due  to  ionospheric  dynamics  and  are  not 
a  physical  characteristic  of  the  ionopause. 


The  outbound  ionopause  crossing  is  poorly  defined 
because  of  a  supertherraal  ion  region  above  270  km  (Taylor 
et  al.,  1980).  As  with  Orbit  170,  a  change  of  the  outbound 
ionopause  altitude  by  about  20  km  will  not  affect  the  anal¬ 
ysis.  At  most,  the  ionopause  crossing  time  will  be  reduced 
by  a  few  seconds,  and  the  only  other  change  would  be  a  small 
decrease  in  the  size  of  the  stable  region  just  below  the 
outbound  ionopause. 
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Figure  3-20.  Orbit  175 — horizontal  projection  of  magnetic 
parameters.  This  is  the  first  of  three  figures  showing 
the  average  field  orientation  every  6  minutes  as  the 
field  rotates  at  a  constant  rate  of  2.5  degrees  per 
minute.  The  valid  altitude  range  is  shown  on  the  right 
hand  edge  of  the  figure  (450-300  km).  The  field  lines 
are  shown  by  the  solid  curved  lines.  The  stability 
contours  (dash  lines)  are:  inner  contour  (shaded  is  at 
350  km  and  the  outer  contour  is  at  450  km.  The  iono¬ 
spheric  stable  region  is  enclosed  by  a  "box."  The 
ionospheric  altitude  range  is  shown  by  the  sidebar;  the 
sidebar  arrows  mark  the  inbound  ionopause  crossing  and 
the  end  of  the  valid  altitude  ranges  (300  km).  The 
observed  magnetic  field  vectors  (arrows)  are  shown  for 
the  ionosheath  and  ionopause  stable  region.  The  cross 
denotes  the  subsolar  point.  The  scale  on  the  upper 
left  hand  corner  of  the  figure  is  the  distance,  in 
degrees,  from  the  plasma  flow  origin  (solid  circle). 
The  crosshatch  area  designataes  the  boundary  of  the 
planetary  disk.  The  asymmetry  factor  is  2.0  and  the 
maximum  field  strength  is  160  gamma. 
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Figure  3-21.  Orbit  175/2 — horizontal  projection.  This 
figure  shows  the  average  field  orientation  6  minutes 
after  Figure  3-20.  The  valid  altitude  range  is  250  km 
down  to  145  km  (periapsis).  The  stability  contours 
are:  inner  contour  (shaded)  is  for  250  km  and  the 

outer  contour  is  for  300  km.  The  sidebar  is  used  to 

i.  .icate  altitudes  within  the  ionosphere.  No  field 
vectors  are  shown  because  this  area  of  the  ionosphere 
is  unstable.  The  plasma  flow  origin  remains  fixed. 
The  asymmetry  factor  is  2.0  and  the  maximum  field 
strength  is  160  gamma. 
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Figure  3-22.  Orbit  175/3 — horizontal  projection.  This 
figure  shows  the  average  field  orientation  6  minutes 
after  Figure  3-21.  The  valid  altitude  range  is  161  km 
up  to  450  km.  The  stability  contours  are:  inner 

contour  (shaded)  is  at  250  km  and  the  outer  contour  is 
at  350  km.  The  ionospheric  stable  regins  are  enclosed 
by  the  shaded  "boxes";  the  shading  indicates  that  the 
stable  region  is  encountered  near  the  250  km  level. 
The  outbound  ionopause  crossing  is  at  300  km.  The 
asymmetry  factor  is  2.0  and  the  maximum  field  strength 
is  160  gamma. 
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Figure  3-23.  Orbit  175 — orbital  plane  cross-section.  The 
contours  of  constant  magnetic  field  strength  are  solid 
lines  in  the  ionosphere  and  dash  lines  in  the  iono- 
sheath.  Contour  magnitudes  are  relative  to  the  stag¬ 
nation  point  maximum.  The  cross  marks  on  the  orbit  are 
separated  by  one  minute  intervals  relative  to  peri- 
apsis.  This  figure  is  a  composite  of  Figures  3-20,  3- 
21,  and  3-22;  stable  regions  are  omitted  because  of 
projection  problems  (see  text). 
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Figure  3-24.  Orbit  175 — comparison  of  model  with  obser¬ 
vations.  Model  results  (dash  lines)  are  shown  only  for 
the  stable  regions.  A  break  in  the  dash  line  design¬ 
ates  a  region  of  instability  predicted  by  the  model. 
The  field  strength  scale  (upper  figure)  is  in  gammas. 
The  asterisks  (*)  are  the  OIMS  measurements  of  atomic 
oxygen  ion  densities.  The  arrows  beneath  the  time 
scale  mark  the  ionopause  crossings.  The  bottom  figure 
shows  how  the  magnetic  field  direction  changes  with 
time. 


TIME  RELATIVE  TO  PERIAPSIS,  min 


4. 


DISCUSSION 


4.1  Comparison  of  Orbits 

Orbits  176  and  177  have  nearly  the  same  magnetic  field 
orientation  relative  to  the  orbital  path,  yet  there  is  a 
significant  difference  in  the  PV  magnetic  field  measurements 
for  these  two  orbits  (see  Figures  3-10  and  3-13).  The  dif¬ 
ference  in  the  observed  magnetic  field  strength  is  most 
likely  due  to  an  observed  decrease  in  the  solar  wind  dynamic 
pressure  (Elphic  et  al .  ,  1980a).  Variability  in  the  solar 
wind  dynamic  pressure  is  also  responsible  for  the  change  in 
size  and  shape  of  the  region  of  stability  for  Orbit  177 
compared  to  Orbit  176  (see  Figures  3-3  and  3-11).  That  is, 
increased  solar  wind  dynamic  pressure  enhances  the  plasma 
depletion  mechanism  which,  in  turn,  increases  both  the 
magnetic  field  strength  and  the  velocity  asymmetry.  Since 
the  ionospheric  magnetic  field  is  closely  coupled  with  the 
plasma  flow,  enhanced  plasma  depletion  also  produces  a 
larger  magnetic  field  asymmetry.  The  change  in  the  field 
asymmetry  reduces  the  angle  9  between  v  and  B  for  inter¬ 
mediate  azimuthal  angles,  and  together  with  the  increased 
field  strength,  these  two  stabilizing  effects  outweigh  the 
destabilizing  effect  of  the  enhanced  flow  velocity.  There¬ 
fore,  the  PV  orbiter  has  the  best  chance  of  intersecting  a 
stable  region  during  periods  of  high  solar  wind  dynamic 
pressure  because,  on  these  days,  the  region  of  stability  is 
enlarged  (i.e..  Orbit  176  versus  Orbit  177).  In  fact,  this 


statistical  prediction  has  been  verified  by  PV  observations 
(Luhmann  et  al .  1980,  1981b). 

The  magnetic  field  observations  are  also  sensitive  to 
changes  of  the  field  orientation  relative  to  the  orbital 
path.  For  example,  if  the  field  pattern  were  rotated  10° 
counterclockwise  in  Figure  3-11  (Orbit  177),  then  the 
orbiter  would  remain  within  the  stable  region  throughout  the 
ionospheric  transit.  However,  a  10°  clockwise  rotation 
would  put  most  of  the  orbital  path  into  an  unstable  area. 
Therefore,  small  field  rotations  can  have  a  dramatic  effect 
on  the  magnetic  field  measurements.  Unfortunately,  similar 
"looking"  magnetic  field  measurements  may  result  from  enti¬ 
rely  different  field  geometries,  and  this  may  mislead  an 
observer  into  thinking  that  similar  "looking"  field  measure¬ 
ments  represent  the  same  physical  ionospheric  conditions. 
For  example,  the  field  strength  measurements  for  Orbit  170 
(Figure  3-7)  and  Orbit  189  (Figure  3-19)  appear  very  simi¬ 
lar,  but  comparing  Figures  3-5  and  3-17  shows  that  the 
ionospheric  conditions  for  these  orbits  were  very  different. 

Orbits  177  (Figure  3-11)  and  188  (Figure  3-14)  provide 
another  useful  example.  The  field  geometry  for  these  two 
orbits  is  similar  but  the  measured  field  strengths  are  very 
different  (see  Figures  3-13  and  3-16).  Certainly  one  dif¬ 
ference  is  that  the  field  geometry  for  Orbit  188  is  rotated 
slightly  relative  to  Orbit  177.  However,  the  most  important 
feature  is  that  periapsis  for  Orbit  188  is  relatively  far 
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from  the  field  pattern  center  and  therefore.  Orbit  188  exits 
the  ionosphere  before  it  can  enter  the  region  of  stability. 
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4.2  Comparison  with  the  2-D  Electrodynamic  Model 

The  modified  two-dimensional  electrodynamic  model 
(Cloutier  et  al ,  1981)  assumed  that  the  ionospheric  magnetic 
field  lines  were  nearly  straight  lines  lying  across  the 
dayside  planetary  disk.  The  two-dimensional  model  horizon¬ 
tal  contours  of  constant  magnetic  field  strength  are  shown 
by  the  solid  lines  in  Figure  4-1  and  the  straight  dashed 

lines  are  the  plasma  flow  streamlines.  The  curved  lines  are 
the  horizontal  field  strength  contours  from  the  three-dimen¬ 
sional  flow/field  model:  only  the  0.6  and  0.2  contours  are 
shown.  Even  though  there  is  a  significant  difference  be¬ 
tween  the  two  models,  the  two-dimensional  model  gives 
excellent  field  strength  results  for  Orbits  175,  176  and 

177.  For  example,  the  two-dimensional  model  analysis  for 

Orbit  177  is  shown  in  Figure  4-2.  Returning  to  the  flow/ 

field  analysis  of  Orbit  177  (see  Figures  3-11,  3-12  and  3- 
13),  notice  that  the  orbital  path  is  nearly  parallel  to  the 
magnetic  equator.  For  this  range  of  azimuthal  angles,  the 
flow/field  model's  magnetic  field  strength  contours  are 

nearly  straight  lines.  Therefore,  although  not  completely 
physically  correct,  the  two-dimensional  model  gives  quanti¬ 
tatively  correct  results  because  the  straight  line  contour 
assumption  was,  in  fact,  a  good  approximation  for  Orbit 
177.  A  similar  argument  is  true  for  Orbits  176  and  175. 

Another  feature  of  the  two-dimensional  model  was  that 
the  instability  criterion  was  written  as: 
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where  C  was  an  empirical  factor  which  ranged  in  value  from 
.01  to  1.0.  The  C  parameter  was  needed  to  explain  the  daily 
variability  of  the  Pioneer-Venus  magnetic  field  strength 
observations.  It  turns  out  that  the  C  parameter  is  nothing 
more  than  a  scaling  factor  which  is  needed  to  properly  pro¬ 
ject  the  three  dimensional  stable  region  geometry  onto  the 
two-dimensional  analysis  plane.  Therefore,  with  the  more 
complete  three-dimensional  flow/field  model,  the  C  parameter 
becomes  unnecessary ,  and  the  stability  ratio  is  simply  given 
by  equation  (2.5.2). 
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Figure  4-1.  Comparison  of  horizontal  magnetic  field 
strength  contours.  The  solid  lines  are  the  two-dimen¬ 
sional  (electrodynamic)  model  contours;  the  contours 
are  labeled  relative  to  the  maximum  field  strength  at 
the  stagnation  point  (i.e.  0.5  =  50%  of  Bmax).  The 
straight  dash  lines  are  the  plasma  flowlines,  and  the 
circle  represents  the  "sub-flow"  region.  The  curved 
lines  are  field  strength  contours  from  the  three  dimen¬ 
sional  model  (only  the  0.6  and  0.2  contours  are  shown). 
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Figure  4-2.  Comparison  of  Pioneer  Venus  observations  with 
the  two  dimensional  electrodynamic  model.  The  upper 
panel  contains  the  observations  and  the  lower  panel 
shows  the  model  results.  The  horizontal  scales  are 
universal  time  (UT)  and  time  relative  to  periapsis 
(RT).  The  bracket  beneath  the  time  axis  indicates  the 
extent  of  the  stable  magnetic  field  region  computed  by 
the  model  (from  Cloutier  etal. ,  1981). 


4.3  Problem  Areas 


One  of  the  major  computational  problems  was  to  maintain 
a  reasonable  decrease  in  the  magnetic  field  strength  near 
the  magnetic  equatorial  terminator.  This  problem  was 
overcome  by  introducing  a  weak  radial  magnetic  field  which 
was  confined  to  a  small  region  about  the  magnetic  equator  at 
large  solar  zenith  angles.  This  weak  radial  field  component 
satisfies  all  the  MHD  equations,  including  Gauss's  Law,  and 
it  contributes,  at  most,  5%  of  the  total  field  strength. 
Therefore,  the  horizontal  assumption  is  still  valid.  A  more 
complete  discussion  is  given  in  Appendix  2. 

Another  problem  is  that  the  model  assumes  that  the 
ionization  is  symmetric  about  the  sub-flow  point  rather  than 
the  subsolar  point.  In  the  dayside  ionosphere,  this  assump¬ 
tion  introduces  only  a  small  error  because  the  ionization 
rate  is  not  very  sensitive  to  small  changes  in  the  solar 
zenith  angle  near  the  subsolar  point  (Chamberlain,  1978). 
Moving  the  ionization  center  to  the  subsolar  point  would 
require  a  major  revision  of  the  fj.ow/field  model  because 
improved  chemistry  algorithms  would  be  needed  to  handle  the 
asymmetry  in  the  ion  production  and  loss  process.  Such  a 
modification  might  be  done  at  a  later  date. 


5. 


CONCLUSION 


The  assumption  that,  far  from  Venus  sub-flow  point,  the 
ionospheric  plasma  flow  and  magnetic  fields  are  strictly 
horizontal  made  it  possible  to  find  analytic  solutions  to 
the  three  dimensional  MHD  equations.  It  was  also  shown  that 
an  asymmetry  exists  in  plasma  flow  distribution:  the  plasma 
flows  faster  perpendicular  to  the  magnetic  fields  than 
parallel  to  these  fields.  This  velocity  asymmetry  was 
largest  near  the  ionopause  and  became  less  pronounced  at  low 
altitudes  because  of  frequent  collisions  with  ambient 
thermal  ions.  High  altitude  plasma  flow  velocities  are  on 
the  order  of  1  to  2  km/sec  and  the  velocity  asymmetry 
(vi/V|-)  can  be  as  high  as  3  to  1.  Since  the  ionosphere  is  a 
good  conductor,  the  magnetic  diffusion  time  is  very  long, 
and  the  close  coupling  of  the  magnetic  field  with  the  plasma 
flow  produced  an  asymmetry  in  the  magnetic  field  geometry. 
Unlike  the  plasma,  the  magnetic  field  asymmetry  is  nearly 
constant  with  altitude.  It  was  also  shown  that  a  velocity 
shear  within  the  plasma  flow  makes  the  ionosphere  suscepti¬ 
ble  to  the  Kelvin-Helmholtz  instability.  This  instability 
occurs  whenever  the  plasma  flow  energy  exceeds  the  restrain¬ 
ing  energy  supplied  by  the  magnetic  field  tension.  Calcula¬ 
tions  showed  that  the  flow  energy  and  field  tension  were 
both  strong  functions  of  azimuth  and  the  resulting  region  of 
stability  has  a  rather  unusual  topology.  This  topology 
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played  a  key  role  in  understanding  the  daily  magnetic 
observations  of  Pioneer-Venus. 

Some  important  results  of  this  study  are: 

1)  the  same  model  can  explain  the  total  spectrum  of 
ionospheric  conditions  ranging  from  total  instability  to 
complete  stability  with  only  small  changes  in  input 
parameters ; 

2)  solar  wind  characteristics  (IMF  direction  and  ram 
pressure)  are  important  factors  in  describing  the  observed 
magnetic  field  observations; 

3)  magnetic  field  strength  measurements  which  appear 
similar  may  represent  very  different  physical  conditions; 

4)  conversely,  similar  magnetic  field  geometries  may 
produce  very  different  magnetic  field  observations  because 
of  small  orbital  alignment  differences; 

5)  the  model  predicts  that  during  periods  of  high  solar 
wind  dynamic  pressure,  the  orbiter  has  a  better  chance  of 
encountering  stable  magnetic  fields  within  the  ionosphere; 
and 

6)  strong  magnetic  field  gradients  are  due  to  iono¬ 
spheric  dynamics  and  are  not  a  unique  physical  character¬ 
istic  of  the  ionopause.  That  is,  the  definition  of  the 
ionopause  most  consistent  with  the  physical  mechanisms  con¬ 
trolling  the  ionosphere  of  Venus  is  the  sharp  gradient  in 
thermal  plasma  marking  the  transition  from  planetary  ioniza¬ 
tion  to  post-shock  solar  wind  plasma. 


APPENDIX:  ANALYTIC  EQUATIONS  FOR  THE  THREE  DIMENSIONAL 


FLOW/FIELD  MODEL 


A. 1 .  Magnetic  Field  Equations 

This  section  contains  the  derivation  of  the  magnetic 
field  equations  discussed  in  Section  2.3  (equations  2.3.1a, 
2.31b  and  2.3.2).  As  stated  earlier,  the  third  equation  of 
motion  (2.1.11)  together  with  Gauss's  Law  (2.1.12)  complete¬ 
ly  describes  the  magnetic  field  geometry.  The  method  of 
solution  is  complicated  because  these  coupled  differential 
equations  do  not  separate  easily.  As  a  matter  of  conven¬ 
ience,  the  third  equation  of  motion  and  Gauss's  Law  are 
repeated  below 


dB_ 


B  f  1  m 
nrsmra  d<li  " 


dB, 

dm 


tanm 


<1  _  _  o  dP 
'  ~  sinm  dt|; 


(A. 1. 1) 


dB 

dm 


m 


B. 


m 


dB 


'l'  _ 


tanm  sinm  d<|> 


=  0 


(A. 1.2) 


where 


B 


Bm  m  +  (|/ 


(A. 1.3) 


The  angles  m  and  if,  and  the  unit  vectors  m  and  <t>  are  shown 
in  Figure  2-1. 

The  first  step  is  to  find  the  homogeneous  solution 
to  the  above  set  of  coupled  equations.  Physically,  the 
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homogeneous  solution  corresponds  to  the  form  of  the  field 
equations  in  the  limit  of  small  <|>  (azimuthal  angle).  As 
discussed  in  Section  2.3,  for  small  4 >,  the  magnetic  field  is 
nearly  parallel  to  4*  and  since  the  plasma  is  free  to  move 
parallel  to  the  field,  the  azimuthal  pressure  gradient  (^0 
will  be  minimized.  Combining  equations  (A. 1.1)  and  (A. 1.2), 
in  the  limit  of  no  azimuthal  pressure  gradient,  we  find  that 


sin2m 


32B 


m 


&B 


am'1 


+  3sinm[cosm  -g^p]  -  Bm  [sin2m-cos2m]  = 


-a23 


m 


m 


54»  ‘ 


(A. 1. 4) 


where  equation  (A. 1.4)  is  in  terms  of  only  one  field  compo¬ 
nent,  Bm.  Next,  using  the  separation  of  variables 

Sm(m,  40  =  Mm)  1(40  (A. 1.5) 

equation  (A. 1.4)  becomes 

sinim  aik  +  3si_nmcosm  ak  _  [ sin2m_cos2m]  =  M  ifl 
k  am2  k  am  I  a<i>2 

(A. 1.6) 

Since  the  right  hand  side  is  in  terms  of  1(40  alone,  then 

-  1  111  =  a2 
2  a*2 


(A. 1.7) 


I 


I  =  c  e 


±ia<J) 


(A. 1. 8) 


where  a  and  c  are  constants.  Substituting  (A. 1.8)  back  into 
(A. 1.6)  and  rearranging  terms  we  find 


tsin2m  kl  -  tSHm  §5  (sin2n  k)  +  U  '  a2]  - 


Using  the  following  change  of  variables 


(A. 1.9) 


f  =  sin2m  k(m) 


z  =  cosm 


(A. 1.10a) 


(A. 1. 10b) 


equation  (A. 1.9)  becomes 

32  f 

(l  -  z2)  -  +  -  (l  -  a2)  =  0 

6z2  (l-z2) 

Possible  solutions  for  (A. 1.11)  are  of  the  form 


(A. 1.11) 


f:(z)  =  Z  trx(l  -  zm) 


nnn  +  X. 


(A.l. 12) 


9 


and 
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f2(z)  =  S  zn  +  X  ) A. 1 . 13 ) 

We  will  look  at  each  of  these  possible  solutions  separately 
and  combine  them  later  to  produce  the  general  homogeneous 
solution. 

Substituting  (A. 1.12)  into  (A. 1.11)  and  solving  for 
various  values  of  m,  we  obtain  the  following  constraints: 


m  =  arbitrary, 

a  =  1, 

n  + 

>> 

II 

O 

(A. 1. 14a) 

m  =  1, 

a  =  1, 

n  + 

\  =  1 

(A. 1.14b) 

m  =  2 , 

o 

II 

n  + 

X  =  1/2 

(A. 1.14c) 

Combining  (A. 1.12)  with 

(A. 1.8) 

and 

(A. 1.5), 

we  find  that 

for  constraint  (A. 1.14a) 


1 

Bra  =  - [c.cos*  +  c  sim|) ]  (A. 1.15) 

a  sin^m  1  * 


for  constraint  (A. 1.14b) 


Bm.  =  - =-  (l  -  cosm)  [c  cos*  +  c  sin* ]  (A. 1.16) 

D  sin^m  3  4 


smm 


(A. 1.17) 
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The  general  solution  for  the  first  case  (A. 1.12)  is  simply  a 
linear  combination  of  (A. 1.15),  (A. 1.16)  and  (A. 1.17): 


B. 


m, 


smm 


,bi  +  b2  cosm. 
4>  +  sin(l>  ( - r-= - J 


sin2m 


(A. 1.18) 


where  bQ,  b:  and  b2  represent  a  collection  of  the  constants 

from  (A. 1.15),  (A. 1.16)  and  (A. 1.17).  In  the  derivation  of 

(A. 1.18),  we  have  applied  the  boundary  condition  that 

1  im  B  =0 
m 

1 

0  (A. 1.19) 


The  constraints  for  the  other  possible  series  solution 
(A. 1.13)  are  found  by  substituting  (A. 1.13)  into  (A. 1.12). 
The  resulting  constraints  are 


a  =  1 ,  n  +  \  =  0 


(A. 1. 20a) 


a  =  1,  n  +  \  =  1 


) A. 1 . 20b) 


and  the  solution  for  Bm2  is 


m. 


(b3  +  b4  cosm) 

-  sin<j» 


sin'-m 


(A. 1. 21) 


where  again  we  have  applied  boundary  condition  (A. 1.19). 
Notice  that  the  solution  (A. 1.21)  is  already  contained  in 
(A. 1.18)  and  therefore,  the  general  solution  for  the 
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component  Bm  is 


bo<t> 

Bm  “  sinm 


+  (b.  +  b,  cosm)  (A. 1.22: 

sin^m  1  i 


A  similar  analysis  can  be  done  for  the  component  and 
the  resulting  solution  is 


ao  cos  cl  r  \ 

B.  =  — -T— -  +  (a  +  a.  cosm  (A. 1.23) 

<1)  smm  e  i  r.  2  m  v  1  2  ‘ 


sin‘m 


where  we  have  applied  the  boundary  condition 


lim 

4>  it/  2 


B,  =  0 


(A. 1. 24) 


Equations  (A. 1.22)  and  (A. 1.23)  can  be  simplified  further  by 
remembering  that  these  equations  must  satisfy  both  the  third 
equation  of  motion  and  Gauss's  Law.  This  additional 
constraint  requires  that  the  constants  in  (A. 1.22)  and 
(A. 1.23)  must  be  related  such  that 

3m  =  U,  -  a.  cosm)  (A. 1.25) 

“H  sm^m  1  1 

=  (a,  -  a,  cosm)  (A. 1.26) 

sin2m  1  2 


The  subscript  H  is  used  to  identify  these  as  the  homogeneous 
solutions  which  are  valid  for  small  <|».  In  order  to  extend 
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the  solutions  for  all  4»,  we  must  return  to  the  inhomogeneous 
coupled  equations  (A. 1.1)  and  (A. 1.2). 

In  order  to  simplify  the  analysis,  we  will  assume  that 
the  azimuthal  pressure  gradient  effects  can  be  completely 
contained  within  the  Bm  component.  Justification  for  this 
assumption  is  both  mathematical  and  physical.  First, 
mathematically  the  assumption  allows  one  to  find  an  analytic 
solution  to  this  coupled  set  of  inhomogeneous  differential 
equations.  This  certainly  isn't  reason  enough  to  use  the 
above  assumption.  However,  the  assumption  produces  analytic 
solutions  whose  form  fits  the  behavior  one  would  physically 
expect.  That  is,  the  field  lines  are  elliptically  shaped 
for  small  <v,  and  then  as  4,  increases,  the  field  lines  turn 
and  become  parallel  to  the  external  IMF  direction.  There¬ 
fore,  we  will  assume  that 


B 


m 


+  B 


mt 


(A. 1. 27) 


B 


(A.l. 28) 


where  BmQ  is  the  term  which  parameterizes  the  effects  of  the 
azimuthal  pressure  gradient.  Using  equations  (A. 1.27)  and 
(A. 1.28),  the  third  equation  of  motion  (A. 1.1)  becomes 


Bm  ( 


_1 _ 

sinm 


3B 


mo 


34, 


_ 1_  3P 

sinm  3 (t> 


(A. 1. 29) 
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and  Gauss's  Law  becomes 


3B, 


m. 


m. 


3m 


tanm 


=  0 


(A. 1.30) 


We  will  further  assume  that  is  separable  and  of  the  form 

Bm  <l>)  =  x  (m)  y  ( <1 )  (A. 1.31) 

0 

Substituting  (A. 1.31)  into  (A. 1.30)  and  integrating,  one 
f  inds 


x(m) 


sinm 


(A. 1.32) 


where  aQ  is  an  integration  constant.  The  form  of  y(<M  is 
more  difficult  to  determine.  First  we  need  to  examine  the 
behavior  of  the  pressure  gradient  term.  As  stated  earlier, 
the  pressure  gradient  goes  to  zero  in  the  limit  of  small  <)/ 
because  the  pressure  is  nearly  uniform  for  small  4».  From 
equation  (A. 1.29)  we  can  see  that  Bmo  itself  must  go  to  zero 
for  small  <|>.  A  form  of  y(4>)  which  satisfies  these 
constraints  is 

ae 

y<4>)  =  Z  c  (<1>  -  sin(|/)Y  (A. 1.33) 

Y=1  Y 

We  will  only  use  the  first  term  of  this  series  because,  even 
with  only  the  first  term,  the  strength  of  the  pressure  grad- 
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ient  is  occasionally  overestimated.  Therefore,  combining 
(A. 1.33)  and  (A. 1.32)  with  (A. 1.31),  BmQ  becomes 


mQ  smm 


(<|i  -  sin<|») 


(A.l. 34) 


Bm  =  sThm  “  sin^  +  ^a2  "  a !  c°sm )  (A.l. 35) 


B0,  =  ~  2"  ^ai  “  a  cosm)  (A. 1.36) 

f  sm^m  1  *■ 


where  the  constants  aQ,  a;  and  a2  are  evaluated  in  Appendix 
A.  2.  This  evaluation  process  is  quite  involved  because  all 
the  necessary  boundary  conditions  are  not  known  and  there¬ 
fore,  we  must  rely  on  geometric  and  physical  arguments  to 
evaluate  the  constants. 

In  practice,  the  field  line  geometry  produced  by 
equations  (A. 1.35)  and  (A. 1.36)  is  correct.  Unfortunately, 
the  field  strength  falls  off  too  slowly  along  the  magnetic 
equator.  That  is,  we  would  expect  that  the  field  strength 
should  decrease  near  the  terminator  at  a  rate  required  by 
the  external  solar  wind  dynamic  pressure.  This  was  verified 
by  Vaisberg  et  al .  (1980a,  1980b)  who  found  that  the  ob¬ 
served  field  strength  variation  was  proportional  to  cos2x» 
where  x  is  the  expected  angle  between  the  normal  to  the 
ionopause  and  the  solar  wind  direction.  As  explained  in 
Section  2.3,  eventhough  Vaisberg  et  al .  based  their  con¬ 
clusion  on  unnormalized  data,  their  proportionally  relation 


135 


is  probably  quite  reasonable.  The  model's  field  strength 
calculations  approximate  the  cos2x  behavior  with  the  intro¬ 
duction  of  a  very  weak  radial  magnetic  field  component. 
This  radial  component  is  confined  to  the  vicinity  of  the 
magnetic  equator  and  rarely  contributes  more  than  5%  of  the 
total  field  strength.  Therefore,  the  horizontal  field 
approximation  is  still  valid.  The  form  of  the  radial 
component  is 

ao,  , 

Br  *  TEamn  t*  “  s1"*)  pm  Vn>  U.1.37) 


where 


P  =  (  1  +  3sin2m)  e~1/(2sin2m) 

m  v  ' 


(A. 1.38) 


P„.(n)  =  sin%  e-1/(sin^)! 


(A. 1. 39) 


The  introduction 

of 

this  weak  radial 

component 

requires  a 

modification  of 

the 

3m  component 

so 

that 

Gauss's  Law  is 

still  obeyed. 

The 

final  form 

of 

the 

field 

equations 

(including  A. 1.37)  becomes 
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A. 2  Determination  of  Constants 


The  first  step  is  to  determine  an  equation  which 
relates  the  total  field  strength,  B,  to  the  angle  0,  where  0 
is  the  angle  between  v  and  B.  As  shown  in  Section  2.3,  the 
magnetic  field  can  be  written  as 

B  =  Bm  m  +  <1  (A. 2.1) 


where 


B_  =  B  cos9 
m 


(A. 2. 2) 


3^  =  B  sin9  (A. 2. 3) 

Substituting  (A. 2. 2)  and  (A. 2. 3)  into  the  third  equation  of 
motion  (A. 1.1)  yields 
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and  Gauss's  law  becomes 
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Combining  (A. 2. 4)  with  (A. 2. 5)  gives 

_1 _  3_  _  B  39  3P 

sinm  34.  2  3m  smm  5cjJ 


(A. 2.6) 


amd 


3_  §2  B2  B2  39  =  * otanQ  3P 

3m  2  tanm  smm  34.  sinm  34.  '  *  * 

Equations  (A. 2. 6)  and  (A. 2.7)  relate  the  magnetic  field 
strength  (B)  to  the  field  line  geometry  specified  by  the 
angle  9.  A  useful  integratable  equation  is  obtained  in  the 
limit  of  small  4»,  where 


lim 

4.  *  0 


3  P 
34> 


0 


(A. 2. 8) 


Therefore,  equations  (A. 2. 6)  and  (A. 2. 7)  become 
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AnB 
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(A. 2.9) 


and 


3 

5m 


InB  + 


JL 3  _  _1 39 

tanm  ~  -  smm  57J7 


(A. 2. 10) 


In  the  limit  of  small  4>,  the  field  line  shape  can  be 
approximated  by  an  ellipse  (see  Figure  A-l).  One  can  show 
that  for  the  spherical  triangle  formed  by  the  tangent  line, 
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the  arc  length  dQ  and  the  arc  length  m'. 


tanQ 


(tan<J>  cosm 


tan((do-d)/R)^_1 
sin4*  sinm  ' 


(A. 2. 11) 


where  the  arc  length  d  is  shown  in  Figure  A-l,  and  R  is  the 
radial  distance  (from  center  of  planet).  Equation  (A. 2. 11) 
is  only  valid  in  the  limit  of  small  4>.  The  second  term  on 
the  right  hand  side  of  equation  (A. 2. 11)  needs  further  eva¬ 
luation  because,  in  the  limit  of  small  (|>,  d  approaches  dQ. 
It's  easiest  to  continue  this  analysis  in  the  projection 
plane . 

In  the  projection  plane 


0o  -  D 


r 


a2  sin24; 

b2  COStJ; 


(A. 2. 12) 


where  r  is  the  radial  distance  (in  the  plane),  a  is  the 
semi-major  axis  and  b  is  the  semi-minor  axis  of  the  ellipse 
in  the  projection  plane,  DQ  is  the  projection  of  dQ,  and  D 
is  the  projection  of  d.  The  angle  ^  does  not  change  in  the 
projection.  Equation  (A. 2. 12)  is  derived  for  the  point 
where  the  tangent  line  touches  the  ellipse.  Figure  A-2 
shows  the  relationship  between  the  projection  and  surface 
parameters  for  small  <|>.  One  can  show  that 

d-d  (d-d) 

sin( — a—)  = - R — 


cosm 


(A. 2. 13) 


Figure  A-l.  Surface  coordinates  for  field  line  geometry 
calculations. 
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Figure  A-2.  Cross-section  near  the  magnetic  equator  showing 
the  relation  between  the  surface  variables  and  the 
projection  plane  variables. 
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and  since  d  approaches  dQ  as  ip  goes  to  zero,  then 
becomes 


dQ  -  d  =  (dq  -  d)  cosm 


Substituting  (A. 2. 14)  into  (A. 2. 12)  yields 


a  a  -  a2  ^  sm2tp 

-  d  =  —  r  cosm 

O  .  2  COSip 


where 


r  =  R  sinm 


In  the  limit  of  small  <p , 


tan 


Therefore  using  (A. 2. 17),  (A. 2. 16)  and  (A. 2. 15), 

(A. 2. 11)  becomes 


and 


tan9  = 


tanip  cosmfl  -  — ) 
b2 


36 
3  ip 


2 

cosm  (l  -  — ) 
b2 


(A. 2. 13) 

(A. 2. 14) 

(A. 2. 15) 

(A. 2. 16) 

(A. 2. 17) 

equation 

(A. 2. 18) 


(A. 2. 19) 
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Equation  (A.  2. 18)  and  (A.  2. 19)  are  only  valid  for  small  4«. 
Finally  substituting  (A. 2. 19)  into  (A. 2. 10),  we  find 

Sm  *"B  =  -  tSHm  (A.2.20) 

In  terms  of  surface  arc  lengths 

v  t 

a  =  R  sin  (-=-)  »  R  sin(m  )  (A. 2. 21) 

K  a 

V||t 

b  =  R  sin  (-£-)  =  R  sin(mb)  (A. 2. 22) 


and 


f  =  a/b 


(A. 2. 23) 


where  f  is  the  surface  asymmetry  factor.  If  we  assume  that 
for  small  <1>,  that  m  *  ma,  then  (A.2.20)  becomes 


in  (§-)  =  -  i  J™  sin  '2^  dm  (A. 2. 24) 
V  7  mi  sin2(™) 

Equation  (A. 2. 24)  can  be  solved  analytically  for  a  variety 
of  field  asymmetry  factors.  Furthermore,  along  the 
perpendicular  direction  ( 4>  -  0),  equation  (A. 1.41)  gives 


a„  cosm) 

2 

a2  cosm ) 


sin2  mi 
sin2  m 


(A. 2. 25) 
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Using  (A. 2. 25)  with  (A. 2. 24),  the  constants  and  a2  can  be 
solved  for  relative  to  some  initial  value  of  B  (the  2-D 
model  of  Cloutier  et  al .  is  a  useful  guide  in  selecting 
B  ^ ) .  Therefore  the  two  unknowns,  al  and  a2  have  been 

replaced  by  a  single  unknown  B^,  which  is  a  model  input. 
Table  A-l  contains  a  list  of  the  constants  aL  and  a2  com¬ 
puted  for  various  asymmetry  factors  and  initial  values  of  B. 

In  order  to  solve  for  the  constants  aQ  and  aQ  ,  one 

l 

needs  to  determine  how  the  field  strength  varies  with  solar 
zenith  angle  along  the  parallel  direction  ( *  ir/2). 
Vaisberg  et  al .  (  1980a,  1980b)  report  that  the  best  fit  to 

the  SZA  variation  of  the  observed  ionopause  field  strength 
is 


B( SZA) 

5 - 

o 


COS2X 


(A. 2. 26) 


where  BQ  is  the  nose  value  for  the  field  strength  and  x  is 
the  angle  between  the  ionopause  and  the  solar  wind  direction 
at  a  given  SZA.  However,  (A. 2. 26)  cannot  be  used  directly 
because  the  PV  observations  were  taken  over  a  period  of 
weeks.  During  this  time,  the  PV  orbit  precessed  toward  the 
terminator,  and  there  is  no  way  to  determine  the  location  of 
the  subflow  point  and  the  orientation  of  the  IMF.  There¬ 
fore,  (A. 2. 26)  is  of  limited  value.  The  Vaisberg  et  al . 
study  does  show  one  important  fact:  the  field  strength  near 
the  terminator  can  be  very  high.  Equation  (A. 2. 26)  will 
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prove  to  be  a  useful  guide  in  evaluating  the  constants 

a 0  and  a^  . 

Unlike  the  constants  a!  and  a2,  there  is  no  analytic 

expression  available  to  evaluate  a0  and  a0  directly. 

1 

Instead,  an  iterative  technique  is  used;  it  begins  by 

picking  a  value  for  the  field  strength  near  the  terminator 

and  at  m  =  40°.  This  is  not  a  totally  arbitrary 

procedure.  That  is,  we  know  the  value  of  3  at  the 

point  (<l>  =  0,  ra  =  40)  and  we  require  that 


BU  =  90,  m  =  40)  v  i 
9 COi  =  0,  m  =  40)  >  i 


(A. 2. 27) 


Physically,  equation  (A.  2. 27 )'  requires  that  the  flow  carries 

the  field  along  faster  in  the  perpendicular  direction  than 

in  the  parallel  direction.  As  a  first  guess,  the  above 

ratio  is  set  equal  to  the  asymmetry  ratio  and  the  constant 

a0  is  set  to  zero.  Next  the  field  strength  is  computed  at 
1 

the  terminator  and  is  then  compared  to  the  Vaisberg  et  al . 

results.  If  there  is  a  substantial  difference,  then  a0  is 

1 

used  to  correct  the  field  values.  The  procedure  is  repeated 
for  different  ratio  values  for  (A. 2. 27).  The  correctness  of 
each  iteration  is  based  on  the  computed  field  line  shape. 
As  discussed  earlier,  the  constants  aQ  and  aQ  represent  the 
effects  of  density  gradients  on  the  field  line  geometry.  By 
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selecting  various  combinations  of  a0  and  a0  ,  we  are  effect¬ 
ively  changing  ionospheric  conditions.  Table  A-l  also 

contains  a  set  of  values  for  a0  and  a0  . 

1 


i 

j 

t 

I 

I 

I 


I 

I 


I 


TABLE  A-l.  CONSTANTS 
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A. 3  Contours  of  Constant  Magnetic  Field  Strength 

The  vertical  field  strength  contours  in  the  magnetic 
prime  meridian  plane  were  developed  in  an  earlier  series  of 
papers  (Daniell  and  Cloutier,  1977;  Cloutier  and  Daniell, 
1979).  The  average  field  strength  distribution  was 

calculated  by  a  variational  technique  which  required  that 
the  ionospheric  volume  Joule  heating  be  a  minimum  for  the 
induced  current  configuration.  This  calculation  showed  that 
the  field  strength  contours  were  nearly  straight  vertical 
lines  above  200  km  up  to  the  ionopause.  The  contours  remain 
vertical  even  with  large  changes  in  the  external  IMF 
magnetic  field  strength  (Cloutier  et  al ,  1981).  Along  the 
magnetic  equator,  the  field  strength  contour  shape  is  found 
using  the  generalized  Ohm's  Law. 

In  MRS  units,  generalized  Ohm's  Law  (Spitzer,  1962, 
Krall  and  Trivelpiece,  1973)  is: 

curl  B  .  3  x  curlB 

curl [  (v  x  B)  -  (VP  +  =— - -)]  =  0 

*0  Me  e  ^o 

(A. 3.1) 

where  v  is  the  plasma  velocity,  B  is  the  magnetic  field 
vector,  a-1  is  the  resistivity  tensor,  Pe  is  the  electron 
pressure,  qe  is  the  electron  charge,  ne  is  the  electron 
density  and  is  the  magnetic  permeability  of  free  space. 
Applying  the  condition  that  B  and  v  are  horizontal  requires 


149 


that 


B 


Bjn  m  + 


4< 


(A. 3. 2a) 


v  =  v  m 


(A. 3. 2b) 


where  m  and  4>  are  the  surface  unit  vectors  defined  in  Figure 
2-1.  If  we  also  restrict  the  analysis  to  the  magnetic  equa¬ 
torial  plane  (along  parallel  axis,  =  0)  and  substituting 
(A. 3. 2)  into  (A. 3.1)  we  find 


m0k 


_  C  vn  x  vt  1  =  5—  ("  -1  5 p .  +  LJk. -  - i —  

qene  ^  e  e  |i0  1 R  5m  ^  R  tanm  R  smm 


bF 


m 


h  m  r  3  p  _F^J_  1  5FR] 

+  JZ  L  F<|,  +  Rsinrn 


+  i_  r-  1  3Fri 

V>o  ^  'R  ■5m“J 


(A. 3. 3) 


where 


B  x  curl  B 


P  = 

(A. 

3.4) 

c  — 

We  have  also 

assumed 

that  at  high 

altitudes, 

the 

conductivity 

becomes  very 

large 

( approaches 

' frozen 

in' 

conditions) 

and 

therefore 

,  the 

second  term 

in  equation 

(A. 3. 1)  is 

ignorable . 

The 

electron 

pressure 

was 

approximated 

by 

the  ideal 

gas  law 

,  where  k 

is 

the  Boltzman 

constant. 


The  first  term  in  equation  (A. 3.1)  vanishes 


because  along  the  magnetic  equator  v  is  parallel  to  B.  One 
would  also  expect  that  the  azimuthal  variation  of  the 

electron  density  and  temperature  near  the  magnetic  equator 

an  aT 

0  0 

would  be  small  (i.e.,  a  =»  0)  and  therefore,  equation 

(A. 3. 3)  becomes 


mek 

%ne 


( 


ane  aTe  3ne 
5r  5m  am 


R  rl  3  F4>  ,  f4»  1  aFnu 

iiQ  31  R  tan'm  R  sinm  5<i  J 


+ 


m 

^o 


c  a  F<j> 

L  m 


F4>  ,  1  aFRi 

R  sinm  8<J»  J 


+ 


4» 

^o 


[- 


i  aFR-, 

T*  -5nTJ 


(A. 3. 5) 


In  terms  of  known  solutions  for  B  (A. 1.35  and  A. 1.36),  F  can 
be  written  as 

lira  F  = 


R_ 

^o 


(S-P9-) 


q  2  *  B  3 

jl.)  +  ± _ OLl 

R  R 


sin2m 


(A. 3. 6) 


Therefore,  the  m  component  of  equation  (A. 3. 5)  must  be  zero, 


1 


dFR 


R  sinm  34; 


0 


(A. 3. 7) 


and  by  substituting  (A. 3. 6)  into  (A. 3.7),  one  can  show  that 
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(A. 3. 8) 


where  Bm^.  is  the  field  strength  value  at  the  height 
(measured  from  center  of  planet).  Since  R  and  R^  can  be 
written  as 


R  =  Rv  +  h  and  h  <<  R 


“  Rv  +  ^k 


hk  <<:  Rk 


(A. 3.9) 


where  Rv  is  the  radius  of  Venus;  then  (A.  3. 8)"  becomes 


8 

S 


m 


m 


k 


r1  +  h/Rv  ^ 
+  \/Rv 


=■  1 


(A. 3. 10) 


Therefore,  the  field  strength  does  not  change  significantly 
with  height  and,  in  the  magnetic  equatorial  plane,  the  field 
strength  contours  are  also  vertical  lines.  It  is  not  pos¬ 
sible  to  solve  for  the  field  strength  contours  at  any  other 
location  without  first  solving  the  full  three  dimensional 
MHD  equations.  Since  the  field  strength  contours  are 
vertical  lines  at  the  two  extremes  (4>  *  0  and  4>  =  it/2),  we 
will  make  the  simplifying  assumption  that  the  field  strength 
contours  are  vertical  for  all  azimuthal  angles. 
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